UNCLASSIFIED 


AD  NUMBER 


AD802503 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors; 
Administrative/Operational  use;  Sep  1966. 
Other  requests  shall  be  referred  to  Army 
Tank-Automotive  Center,  Warren  MI. 


AUTHORITY 


USATAC  ltr,  16  Apr  1969 


THIS  PAGE  IS  UNCLASSIFIED 


802503 


TECHNIC 


AL  REPORT  NO. 


9387  Cll  109) 


ATLAS  OF  0FF-R0A0  GROUND  ROUGHNESS  P.S.O.'s  AND 
REPORT  ON  DATA  ACQUISITION  TECHNIQUE 


John  L.  Bogdanoff 
Frank  Kozin 
Louis  J.  Cote 


DEVELOPMENT  LABORATORIES 


U  S.  ARMY  TANK  AUTOMOTIVE  CENTER 

\  .  '  \ 


Tn*  findings  in  this  f  sport  or*  not  to  0* 
conitrusd  os  «n  official  Ospartmsnt  of  in* 
Army  position,  unless  so  <Msignat«4  by  o<n*r 
■uthorusd  document*. 


DDC  AVAILABILITY 


QUALIFIED  REQUESTERS  MAY  OBTAIS 
COPIES  OF  THIS  REPOST  PROM  DDC. 


CITATI0K  OP  EQUIPMENT  Iff  THIS  REPORT 
DOES  MOT  CONSTITUTE  AH  OFFICIAL 
IHDORSEMEHT  OR  APPROVAL  OP  THE  USB  Of 
SUCH  COMMERCIAL  HARDWARE., 


Dwtroy  thh  report  *t*m  ft  to  no  longer 
HmM.  Do  not  return  it  to  th*  origmetor. 


AtfCHS  Codet  5022.11.022  Project  Order  No. i  1-V-0-21701-A-045 

PRON*  A 1 -7-R002 1 -A 1 -£H 

I 

I 

LANO  LOCOMOTION  LABORATORY 


ABSTRACT 


Power  spectral  densities.  In  one (line) and  two (area) 
dimensions,  are  presented  for  off-road  ground  In  eleven 
sites  in  the  United  States.  Methods  of  acquiring,  record 
ing,  and  processing  the  data  are  described  In  detail. 


FOREWORD 


Ground  roughness  measurements  on  lines  at  three 
sites  were  given  in  a  previous  report  and  their  power 
spectral  densities  (p.s.d.'s)  displayed.  These  p.s.d.'s 
possessed  many  common  characteristics  and  the  question 
was  raised  by  many  as  to  whether  these  characteristics 
would  persist  for  other  sites.  Moreover,  no  computer 
program  was  available  for  processing  area  data  or  even 
parallel  line  data. 

These  basic  data  and  their  p.s.d.'s  are  required  in 
any  vibrational  analysis  of  the  motion  of  a  vehicle  on 
open  ground. 

The  opportunity  arose  in  the  summer  of  1964  to 
acquire  additional  roughness  data  in  the  mid-continent 
of  the  United  States. 
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GLOSSARY 


Ground  Profile  Is  a  plot  of  the  survey  height  vs  . 
distance  along  a  line . 

Power  Spectral  Density  (p.s.d.)  measures  the  amount 
of  vibration,  by.  frequency  bands,  of  the  ground  heights. 
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INTRODUCTION 


The  theory  of  land  locomotion  tre  are  developing 
requires  ground  roughness  spectra  as  Input  data.  Line 
spectra  have  been  obtained .  from  a  limited  number  of  sites. 
However,  vehicles  make  use  of  parallel  tracks  and  the 
dynamical  Input  from  parallel  tracks  requires  the  spectra 
on  each  track  and  also  the  cospectra.  A  complete  des¬ 
cription  of  a  site  would  require  the  spectra  of  a  track 
and  the  cospectra  of  parallel  tracks  with  different 
separations.  Also,  the  spectra  and  co -spectra. for  tracks 
in  different  directions  would  be  required.  Fortunately, 
this  total  Information  is  contained  In  the  two-dimensional 
spectrum  [See  Report  [1]  (references)  for  a  complete  des¬ 
cription]. 

The  preparation  of  Inputs  to  the  differential  equa¬ 
tions  of  motion  making  use  of  the  two-dimensional  spectrum 
is  described  in  References  [2].  We  are  now  in  the 
process  of  comparing  this  first  order  theory  with  experi¬ 
ment  . 

Part  of  the  basic  information  needed  In  this  theory 
is  a  description  of  ground  roughness  in  terms  of  p.s.d.'s 
(power  spectral  densities).  In  one  of  our  reports  [1], 
we  presented  p.s.d.'s  obtained  from  data  taken  at  three 
sites.  Two  of  the  sites  had  been  modified  by  vehicle 
traffic.  The  three  spectra  had  many  points  of  similarity 
and  some  points  of  difference.  We  interpreted  the  spectra 
and  compared  them  with  the  profile  data.  Certain  questions 
remained  outstanding  which  required,  for  their  resolution, 
the  acquisition  of  additional  data: 

a)  Are  these  p.s.d.'s  typical? 

b)  Can  the  major  and  relevant  features  of  ground 
p.s.d.'s  be  summarized  In  a  simple  manner? 

c)  Is  It  possible  to  estimate  the  relevant  charac¬ 
teristics  of  ground  roughness  p.s.d.'s  from  infor¬ 
mation  found  In  aerial  photography,  topographical 
and  geographical  maps,  etc.? 

It  was  therefore  decided  to  compile  a  small  atlas  of 
ground  roughness  p.s.d.'s  made  at  a  variety  of  sites. 
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During  a  visit  (May  1964)  to  Vicksburg  the  objec¬ 
tives  were  discussed.  It  was  suggested  there  that  we 
choose  our  sites  at  military  installations  that  were  being 
used  by  the  MERS  Project.  This  had  the  advantage  that 
the  sites  were  available  for  future  vehicle  testing.  The 
MERS  Project  had  compiled  photographs  and  topographical 
maps  of  these  installations,  and  the  commanding  officers 
had  indicated  their  willingness  to  cooperate  with  MERS. 

The  location  and  terrain  characteristics  of  the  instal¬ 
lations  were  reviewed  and  preliminary  selection  was  made 
at  that  time. 

These  selected  installations  were  in  the  mid-continent  • 
so  that  a  surveying  crew  could  travel  from  one  to  another 
in  a  short  summer  period  and  make  a  maximum  number  of. 
surveys  at  a  variety  of  locations. 

The  Commanding  Officer  at  ATAC  sent  a  letter  to  each 
installation  explaining  our  purpose  and  requesting  per¬ 
mission  to  make  an  initial  visit  and  then  to  make  the 
required  measurements  at  a  subsequent  time.  It  was  also 
explained  that  we  might  need  the  services  of  a  guide,  etc. 

Preliminary  trips  were  made  by  MASC  personnel  to 
select  sites  at: 

1.  Port  Riley,  Kansas 

2.  Camp  Gruber,  Oklahoma  (abandoned) 

3.  Fort  Carson,  Colorado 

4.  Port  Sill,  Oklahoma 

5.  Port  Hood  ,  Texas 

6.  Port  Polk,  Louisiana 

7.  Port  Benning,  Georgia 

8.  Port  McClellan,  Alabama 

Two  sites  were  selected  at  each  installation  aft'er  a 
reconnaissance  of  the  general  areas  available  to  us.  The 
sites  were  selected  so  there  would  be  a  maximum  chance  of 
availability  within  the  next  few  years  if  vehicle  tests 
were  desired.  In  addition,  the  following  considerations 
were  used  in  site  selection: 
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a)  relatively  free  of  heavy  ground  cover  so  that 
the  surveying  operation  would  not  be  unduly 
hampered, 

b)  relatively  near  access  roads, 

c)  relatively  flat  and  of  fairly  uniform  rough¬ 
ness  characteristics  over  fairly  large  areas, 

d)  relatively  undeformed  by  vehicles. 


Under  no  circumstances  was  it  very  difficult  to  meet 
the  above  conditions. 

In  all  cases,  the  commanding  officers  and  personnel 
of  the  installations  gave  outstanding  cooperation. 

In  addition  to  the  data  collected  above,  other  survey 
dataare  available  at  MASC.  These  data  will  also  be 
Incl'ded  in  the  atlas.  The  locations  of  these  additional 
sit  -  are  as  follows: 

9.  Port  Knox,  Kentucky 

10.  Aberdeen  Proving  Grounds,  Maryland 

11.  Yuma  Test  Station,  Arizona 

12.  Las  Vegas,  Nevada  (Mercury  Test  Site  of  AEC) 

Vehicle  tests  were  conducted  at  Port  Sill,  Las  Vegas, 
and  Thailand  along  with  elevation  surveys.  Only  the 
survey  data  associated  with  the  first  two  vehicle  test 
sites  will  be  reported  upon  here;  the  vibration  data  will 
be  discussed  in  a  subsequent  report.  The  Thailand  data 
a*.e  on  a  road  rather  than  open  ground,  and  thus  will  not 
appear  in  this  atlas. 

With  respect  to  the  surveying  program,  a  number  of 
decisions  were  made  by  MASC. 

It  was  decided  to  survey  two  sites  at  each  instal¬ 
lation,  preferably  of  different  character.  Each  site 
would  provide  area  and  line  data.  Preliminary  recon¬ 
naissance  Indicated  (3ee  [1]  for  a  discussion  of  details 
in  planning  experiments  of  this  type)  that  a  two-foot 
spacing  would  be  sufficiently  close  to  include  all  sig¬ 
nificantly  high  frequencies  without  aliasing.  (Observa¬ 
tion  of  the  vehicle  trials  at  Fort  Sill  suggests,  however, 
that  closer  spacing  may  be  necessary  on  rocky  ground.) 


The  number  of  data  points  In  a  square  was  decided  to  be 
10,000  (100  x  100)  and  for  the  line  500.  The  dimensions 
of  the  squares,  therefore,  were  200  ft.  x  200  ft.  and  the 
lines  were  1000'  long.  Por  convenience,  we  required  the 
line  to  go  through  the  center  of  the  square.  However, 
the  angle  of  the  line  to  the  side  of  the  square  was  left 
to  the  surveyors  and  a  variety  of  angles  was  expected 
and  obtained  from  the  different  sites . 

Preliminary  discussion  with  John  Chen,  a  professional 
land  surveyor,  indicated  that  normal  horizontal  and  ver¬ 
tical  control  within  1/100  of  a  foot  and  horizontal  con¬ 
trol  2/100  of  a  foot  could  easily  be  obtained.  This  is 
well  within  limits  required  by  the  dynamic  problem.  In 
view  of  the  fact  that  loose  material  on  the  surface  of 
the  ground  will  give  comparable  deviation,  the  control 
was  deemed  adequate. 

A  survey  crew  was  organized  under  the  leadership  of 
Dr.  Jay  Barton,  a  biophysicist  at  St.  Joseph  College, 
Rensselaer,  Indiana.  .  The  crew  consisted  of  Arthur  Hawkin3 
Donald  Paarlberg,  Jr.,  John  Poster,  Paul  F.  Chenea,  Jr.,, 
and  William  Whistler.  Of  these  Dr.  Barton  and  Arthur 
Hawkins  had  surveying  experience.  The  others  were  college 
students  on  summer  vacation. 

Under  the  direction  of  Dr.  Barton,  and  with  advice 
of  Mr.  John  Chen,  equipment,  and  other  necessary  items 
were  selected  and  acquired.'  The  equipment  and  crew  were 
fitted  into  a  standard  Chevrolet  Station  Wagon  rented 
for  the  purpose.  There  was  a  short  training  period 
during  which  the  crew  practiced  stowing  the  equipment 
and  carrying  out  the  surveys.  The  tour  was  completed 
ahead  of  schedule  without  mishaps  to  personnel  or  equip¬ 
ment  . 

On  the  average,  it  took  four  days  to  complete  the 
surveys  at  an  installation  and  approximately  one  and  a 
half  days  to  travel  between  installations.  Approximately 
180,000  data  points  were  acquired  by  this  crew. 

The  surveying  equipment  included  the  following  major 
Items: 

2  Repeating  Theodolites,  K&E  730050 
2  Self -Leveling  Zeiss  Levels,  K&E  750020 
2  Lenker  Elevation  Leveling  Rods 
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In  addition ,  the  equipment  Included  both  steel  and  cloth 
tapes,  marking  pins,  range  poles,  etc.  A  complete  list  of 
surveying  equipment  is  given  in  Appendix  A. 


In  addition  to  the  survey  data,  the  following  data 
were  acquired  at  each  site.  The  reconnaissance  party  lo¬ 
cated  the  survey  sites  on  topographical  maps,  and  marked 
them  on  the  ground' wlthr  paint,  plastic  tape,  etc.  The 
orientation  of  the  survey  lines  and  squares  was  noted  by 
the  survey  party.  Verbal  descriptions  were  made  and  in 
some  cases  photographs  were  taken  to  further  Identify 
sites. 


The  survey  data  plus  the  other  material  discussed 
above  taken  at  a  site  were  returned  each  week  to  MASC 
by  registered  mail.  On  his  return  from  the  surveying 
operations,  Dr.  Barton  put  the  survey  data  in  order  for 
data  processing.  Dr.  Barton  also  remained  on  the  staff 
of  MASC  to  assist  in  the  data  processing. 


The  data  were  put  on  IBM  cards ,  the  cards  were  pre- 
p  ocessed  to  remove  the  Instrument  height  and , compare  each 
aata  point  with  Its  nearest  neighbors.  The  output  was 
listed  and  large  deviations  were  marked  by  an  asterisk. 
These  marked  points  were  compared  by  Dr.  Barton  with  the 
data  in  the  survey  books  and  any  errors  were  corrected. 

The  corrections  were  not  made  on  the  IBM  cards  but  rather 
on  the  input  tapes.  Both  cards  and  tape  are  stored  at 
MASC  and  copies  are  available  at  cost  upon  request. 

This  report  is  the  first  of  two  reports.  In  this 
report,  we  shall  be  concerned  with  presenting  the  com¬ 
putational  results,  the  site  description,  and  the  com¬ 
puting  programs.  Conclusions  and  recommendations  will 
be  concerned  with  data  acquisition  methods  and  process¬ 
ing.  The  second  report  will  deal  with  the  interpreta¬ 
tion  and  implication  of  these  results. 


OBJECT 


The  object  of  this  report  is  to  present  the  results 
of  mid-continent  surveys  and  the  computations  performed. 
Interpretations  will,  to  a  large  part,  be  presented  in  the 
next  report.  Recommendations  will  be  limited  to  the  record 
ing,  acquisition,  and  processing  of  the  data. 


SUMMARY 


A  short  description  is  given  of  the  survey  methods  and 
the  equipment  used.  This  is  followed  by  a  table  listing 
all  site  locations  and  orientations  of  the  squares  and  line 
surveyed.  The  results  (p.s.d.’s  in  numerical  form  and  graphs 
of  the  p.s.d.'s)  of  one-dimensional  (line)  surveys  are  pre-y 
sented  next;  results  for  parallel  lines  are  also  included;" 
The  section  which  follows  presents  similar  results  for  two- 
dimensional  (area)  data. 

Complete  descriptions  of  the  two  computer  programs  are 
given. 


CONCLUSIONS 


The  p.s.d.’s  obtained  from  the  data  taken  at  the  addi¬ 
tional  sites  show  the  same  general  features  noted  in  the 
previous  report . 

Rod  and  level  survey  methods  of  data  acquisition  are 
cheap  and  efficient. 
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RECOMMENDATIONS 


Further  data  acquisition  by  the  methods  used  for  the 
purpose  of  enlarging  this  atlas  are  not  recommended. 

When  roughness  data  are  required  for  special  purposes, 
however,  the  rod  and  level  method  used  is  recommended. 
Self-leveling  levels  and  self-zeroing  rods  are  recommended 
as  basic  equipment. 
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SITE  DESCRIPTION 


Survey  sites  are  located  with  reasonable  accuracy  on 
topographical  maps  along  with  verbal  descriptions  ar.d  in 
some  cases  photographs.  No  attempt  was  made  to  preserve 
the  exact  locations  of  data  points;  first,  this  is  not 
relevant  from  the  point  of  view  of  p.s.d.  analysis;  second, 
to  do  this,  permanent  markers  would  have  had  to  be  erected 
and  permission  to  do  this  was  thought  to  be  hard  to  obtain. 

At  this  point,  it  should  be  emphasized  that  the  data  were 

acquired  to  make  a  statistical  estimate  of  a  roughness  char¬ 
acteristic  and  not  to  provide  an  exact  description  of  ground 
elevation.  Our  results  should  be  regarded  as  reproducible 
within  the  statistical  accuracy  to  be  described  later.  For 
this  reason,  the  exact  relocation  of  data  points  is  consid¬ 
ered  to  be  unnecessary. 

If  vehicle  tests  are  to  be  conducted  at  a  site  at  which 

a  survey  has  been  made,  our  results  will  be  useful  in  gen¬ 

eral  selection  considerations.  When  vehicle  tests  are  made, 
our  exact  site  might  not  be  the  most  convenient;  even  if  it 
is,  ground  deformation  may  occur  due  to  a  variety  of  factors 
associated  with  vehicle  weight,  ground  moisture,  etc.  Hence, 
surveys  along  actual  tracks  may  prove  necessary. 

The  site  description  information  is  summarized  in  Table 
I.  Additional  information  is  available  at  MASC. 

The  sites  coded  by  letters  A-U  are  those  surveyed  in  the 
summer  of  1964.  Some  additional  lines  were  surveyed  at  Ft. 

Sill  in  connection  with  vehicle  trials  conducted  there.  Lines 
L  and  M  are  on  half  of  the  previously  surveyed  lines  of  G  and 
H.  In  particular,  care  was  taken  in  conducting  the  trial  at 
Site  2  so  that  the  wheel  track  was  exactly  on  the  1000  foot 
line  and  so  the  first  half  of  line  H  represents  the  ground 
heights  before  and  M  the  ground  heights  after  modification 
by  the  vehicle. 

The  remaining  sites,  given  by  names,  were  surveyed  earlier. 
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LINE  SPECTRA 


Our  methods  of  obtaining  p.s.d.  estimates  for  line 
data  are  described  in  [1].  They  will  not  be  repeated 
here.  However,  a  description  of  computation  is  given 
in  the  section  on  computations. 

We  present  in  this  section  the  spectral  estimates 
in  a  numerical  table,  and  graphs  on  semi-log  paper  of 
the  spectral  estimates.  Each  of  the  two-track  spectra 
will  be  graphed  on  the  same  sheet  for  comparison  pur¬ 
poses. 

Co-spectra  between  parallel  lines  are  needed  in 
vehicle  motion  analysis.  They  have  been  computed  where 
parallel  line  surveys  were  conducted.  These  co-spectra 
are  presented  in  tabular  form  only,  since  plots  are  dif¬ 
ficult  to  interpret. 

The  statistical  accuracy  of  the  spectra  of  the  1000 
feet  lines  may  be  given  in  several  ways.  Each  estimate 
has  18  degrees  of  freedom  or,  In  other  words,  a  confi¬ 
dence  interval  at  95%  confidence  may  be  found  by  using 
the  factors  .571  and  2.186.  [We  are  95%  confident  that 
the  true  spectral  value  lies  between  .571  times  the 
estimate  and  2.186  times  the  estimate.]  The  computations 
are  made  according  to  pp  28-30  of  [1], 

The  raw  elevation  data  on  IBM  cards  may  be  obtained 
at  cost  from  MASC. 


spectral  estimates  for 
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_  _  ........  .  : 

..  . .  - - 

9 

1.836E-02 

K  ,r  ■  ..  ■ 

10 

2.1306-02 

. .  - . 

11 

2.396E-02 

•  ■  •  •'  '.  ,i*v' 

12 

3.677E-02 

„  ...  . 

13 

3.706E-02 

14 

2.054E-02 

.....  . . . 

15 

1. 1476-02 

16 

8.753E-03 

. '■  ■--- 

17 

6.6206-03 

« 

13 

6.884E-03 

• .  ......  - .  -  ... 

. -  —  . . '• 

13 

7.776E-03 

20 

7.851E-03 

. -  -  ■■■■■■. 

21 

7.742Er-03 

•  ' 

22 

8.993E-03 

23 

9.915E-03 

24 

7.609E-03 

....  .  '  . 

. 

25 

5.666E-03 

26 

8.499E-03 

•-  • 

- 

27 

1. 122E-02 

28 

1 .2726-02 

29 

1.267E-02 

30 

1.088E-02 

31 

9.100E-03 

32 

7.5466-03 

.  _  . 

33 

6.229E-03 

34 

4.186E-03 

■  • 

35 

4.337E— 03 

36 

5.787E-03 

37 

4.834E— 03 

38 

4. 550E-03 

.  .  .  -  - 

39 

5.4066-03 

40 

5.327E-03 

41 

5.308E-03 

42 

5.299E-03 

43 

4.685E-03 

44 

4.082E— 03 

45 

4.433E— 03 

46 

6.458E-03 

47 

8.227E-03 

48 

6.7216-03 

49 

5.557E-03 

50 

6.1206-03 

i 

i 


i  < 


j  ; 
\ 


' 

\  • 
» .  • 


RILEY  2  LINE  B 


o 


SPECTRAL  ESTIMATES  FOR 

l 

ITEM  A 


n 

p 

WXCP) 

...  0 

— — 

1 

.  1.203E  00 

....  2 

2.243E-01 

3 

1.087E-01 

4 : 

8.5346-02 

5 

4.978E-02 

6 

2.396E-02 

7 

3.1256-02 

.  8 

5.646E-02 

9 

7.093E-02 

10 

6.9726-02 

11 

6.0496-02 

12 

4.057E-02 

O 

13 

3.073E-02 

14 

4.0876-02 

15 

5.1756-02 

w 

16 

6.3656-02 

17 

6.7716-02 

18 

5.056E-02 

19 

2.8386-02 

_  . 

20 

1.7956-02 

21 

2.265E-02 

— ' 

22 

3.0716-02 

23 

2.921E-02 

24 

2.4386-02 

— ■ 

25 

2.7876-02 

. .  .... 

26 

3.3556-02 

27 

3.0116-02 

28 

2. 4276-02 

. 

29 

2.1036-02 

— .  ■ 

30 

2.4786-02 

31 

3.032E-02 

32 

2.8266-02 

33 

2.8056-02 

_  -  - 

34 

2.6176-02 

35 

2.042E-02 

36 

1.9786-02 

37 

2.1696-02 

38 

1.9986-02 

39 

1.7596-02 

40 

1.6706-02 

41 

1.2696-02 

42 

1.1006-02 

- 

43 

1.791E-02 

44 

2.6036-02 

45 

2. 125E-02 

46 

1.073E-02 

47 

9.8506-03 

48 

1.2006-02 

- 

49 

1. 544E-02 

50 

1.8636-02 

?  .^^F?  .  :T.'  «^v*JMsis{n^-ji 


r-'^-yr  -!"■■■*-  v 


SPECTRAL  ESTIMATES  FOR  GRUBER  1  LINE  C  -33* 

ITEM  A 


p 

rt 

WXCP) 

U 

1 

2.149E  00 

2 

2.8806-01 

3 

6.543E-02 

4 

1.565E-02 

5 

1.711E-02 

j 

2.033E-02 

7 

1.835E-02 

3 

1.607E-02 

9 

1.239E-02 

10 

8.907E-03 

11 

9.608E-03 

12 

1.143E-02 

13 

1.338E-02 

14 

1.534E-02 

15 

1.142E-02 

16 

6. 170E-03 

17 

5.095E-03 

13 

5.554E-03 

19 

6.106E-03 

20 

7.327E-03 

21 

9.803E-03 

22 

9.189E-03 

23 

7.117S-03 

24 

7.438E-03 

25 

6.614E-03 

26 

6.596E-03 

27 

9.409E-03 

28 

.  9.730E-03 

29 

7.815E-03 

30 

8. 589E— 03 

31 

9.228E-03 

32 

6.939E-03 

33 

4.723E-03 

34 

5.850E-03 

35 

8.545E-03 

36 

9.1836-03 

37 

7.529E-03 

33 

6.140E-03 

39 

8.922E-03 

40 

1.0076-02 

41 

8.365E-03 

42 

1.115E-0? 

43 

1.144E-02 

44 

7.1096-03 

45 

4.818E-03 

46 

5.149E-03 

47 

6.473E-03 

43 

6.2936-03 

49 

6.853E-03 

\  i 

\ 

SPECTRAL  ESTIMATES  FOR  _  GRUBER,  2  LINE  0  -34 

_ITEM  A _ ..  . 


p . 

NXCPJ 

c 

— 

1 

1.905E  01 

2 

3.259E  00 

3 

8.8416-01 

4 

3.690E-01 

5 

1.95CE-01 

6 

1.425E-01 

7 

1.285E-01 

8 

I. 1906-01 

9 

1. 185E-01 

10 

9.107E-02 

11 

5.692E-02 

12 

4.2566-02 

13 

4.023E-02 

14 

4.4466-02 

15  . 

4.9986-02 

16 

4.7106-02 

17 

3. 1436-02 

18 

2.917E-02 

19 

3. 7016-02 

2C 

2.718E-02 

21 

1.7566-02 

22 

1.775E-02 

23 

1.6606-0 2 

24 

1.721E— 02 

25 

2.067E-02 

26 

1.911E— 02 

27 

1.374E-02 

28 

1.2556-02 

29 

1.1946-02 

30 

1.0406-02 

31 

9.92CE— 03 

32 

9.679E-03 

33 

1.006E-02 

34 

1.124E-Q2 

35 

1.1636-02 

36 

1.027E-G2 

37 

8.567E-03 

38 

7.415E-03 

39 

6.7756-03 

4C 

6.22CE-03 

41 

7.4336-03 

42 

1.0266-02 

43 

9.497E-03 

44 

6. 0226-03 

45 

5. 512 E- 03 

46 

5.822E-03 

47 

6. 107E-03 

48 

6.282E-03 

49 

4.647E-03 

50 

3.7366-03 

SPECTRAL  ESTIMATES  FOR  CARSON  1  LINE  E 
ITEM  A 


P  MX(P) 


0 

—— 

1 

2.411E  00 

2 

2.971E-01 

3 

1.015E-01 

4 

6.537E-02 

5 

5.168E-02 

6 

3.770E-02 

•  7 

2.962E-02 

8 

2.752E-02 

9 

2.460E-02 

10 

1.871E-02 

11 

1.333E-02 

12 

9.1546-03 

13 

6.0606—03 

14 

7.755E-03 

15 

1.298E-02 

16 

1. 3996-02 

17 

1.167E-02 

18 

1.011E— 02 

19 

8.7646-03 

20 

1.105E-02 

21 

1.597E— 02 

22 

1.445E-02 

23 

7.454E— 03 

24 

5.701E-03 

25 

7.  580E-03 

26 

7.095E-03 

27 

6.790E— 03 

28 

5.732E-03 

29 

3.538E— 03.. 

30 

3.7766-03 

31 

5.034E— 03 

32 

5.237E— 03 

33 

4.777E— 03 

34 

4.040E— 03 

35 

4.127E-03 

36 

4.356E-03 

37 

3.822E-03 

38 

4. 140E-03 

39 

5. 120E-03 

40 

5.476E-03 

41 

5.847E-03 

42 

4.943E-03 

43 

2.794E-03 

44 

2.468E-03 

45 

3.352E-03 

46 

3.665E— 03 

47 

3.471E-03 

48 

3.494E-03 

49 

4.302E-03 

50 

4.853E-03 

*.-35 


SPECTRAL  ESTIMATES  FOR  CARSON  2  LINE  F  -36 

ITEM  A 


p 

WX<P) 

0 

1 

9.254E  00 

2 

8.074E-01 

3 

2.365E-01 

4 

1.093E-01 

5 

5.708E-02 

6 

3.408E-02 

7 

1.997E-02 

a 

1.198E-02 

9 

1.065E-02 

10 

1.071E-02 

11 

1.153E-02 

12 

1.203E-02 

13 

1.184E-02 

14 

1.270E-02 

15 

9.685E-03 

16 

5.107E-03 

17 

5.396E—03 

18 

5.543E-03 

19 

4.689E— 03 

20 

5.954E-03 

21 

6.193E-03 

22 

4.367E-03 

23 

3.330E-03 

24 

4.511E-03 

25 

5.787E-03 

26 

5.838E-03 

27 

5.847E-03 

23 

4.197E-03 

29 

3.455E-03 

30 

4.604E— 03 

31 

3.809E— 03 

32 

2.617E-03 

33 

3.245E-03 

34 

5.835E-03 

35 

7.199E-03 

36 

5.620E-03 

37 

3.893E-03 

33 

3.313E-03 

39 

2.983E-03 

40 

2.369E-03 

41 

1.949E-03 

42 

1.733E-03 

43 

2.866E-03 

44 

4.904E-03 

45 

4.194E-03 

46 

2.886E-03 

47 

4.084E-03 

48 

5.077E-03 

49 

4.192E-03 

50 

3.344E-03 

ai 


i 


SPECTRAL  ESTIMATES  FOR  „  SILL  1  LINE  6 


ITEM  A 


p 

WXCP  J 

0 

■Ml.— 

1 

2.405E  01 

2 

2.761E  00 

3 

1.296E  00 

4 

6.264E-01 

5 

2.741E-01 

6 

1.579E-01 

7 

1.290E-01 

8 

7.999E-02 

9 

3.755E-02 

10 

2.126E-02 

11 

1.836E-02 

12 

2.173E-02 

13 

1.874E-02 

14 

1.608E-02 

15 

1.791E-02 

16 

1.915E-02 

17 

1.858E-02 

13 

1.466E-02 

19 

1.503E-02 

20 

1.496E-02 

21 

1.066E-02 

22 

1.032E-02 

23 

9.479E-03 

24 

8.905E— 03 

25 

1.032E-02 

26 

1.046E-02 

27 

1.253E-02 

28 

1.309E-02 

29 

1.105E-02 

30 

1.301E-02 

31 

1.293E-02 

32 

1.013E-02 

33 

1.229E-02 

34 

1.227E-02 

35 

1. 101E-02 

36 

1.210E-02 

37 

1.250E-02 

38 

1.455E-02 

39 

1.439E— 02 

40 

1.191E-02 

41 

1.003E-02 

42 

8.425E-03 

43 

9.439E-03 

44 

1.121E-02 

45 

1.077E-02 

46 

1.017E-02 

47 

8.973E-03 

48 

6.466E-03 

49 

7.605E-03 

50 

9.717E-03 

SILL  2  LINE  H 


SPECTRAL  ESTIMATES  FOR 
ITEM  A 

P  WX{P>  ' 


0 


1 

3.489E  01 

2 

4.077E  00 

3 

1.342E  00 

4 

7.074E-01 

5 

3.865E-01 

6 

1.853E-01 

7 

1.196E-01 

8 

8.334E-02 

9 

5.525E-02 

10 

5.448E-02 

11 

4.154E-02 

12 

2.404E— 02 

13 

2.607E-02 

14 

2.634E-02 

15 

1.712E-02 

16 

1.632E-02 

17 

1.431E-02 

18 

8.266E-03 

19 

1.049E-02 

20 

1.037E-02 

21 

7.043E-03 

22 

7.252E-03 

23 

6.250E-03 

24 

5.936E-03 

25 

8.731E-03 

26 

8.323E-03 

27 

7.559E-03 

28 

9.244E-03 

29 

6.715E-03 

30 

4.445E-03 

31 

4.353E-03 

32 

3.243E-03 

33 

6.  118E-03 

34 

9.188E-03 

35 

5.943E-03 

36 

4.628E-03 

37 

6.348E-03 

38 

5. 319E-03 

39 

6.480E-03 

40 

8.311E-03 

41 

5.677E-03 

42 

5. 104E-03 

43 

6.289E-03 

44 

4. 532E-03 

45 

4.892E-03 

46 

5.807E-03 

47 

5.060E-03 

48 

6.900E-03 

49 

6.405E-03 

50 

4.085E-03 

r 

- 

ft 


i 


i 


SPECTRAL  ESTIMATES  FOR 
_ ITEM  A  KITH  ITEM  8 


SILL  3  LINE  J  2  TRACK '  -39 


WXIPJ 


WYIPI 


HCCP) 


WQfPI 


'6.9526-01 

6.7916-02 

3.9806-02 

3.3006-02 

1.926E-02 


6 

9.2516-03 

‘7 

7.3376-03  " 

8 

9.8906-03 

9 

1.237E-02 

10 

1.3206-02 

11  ”* 

”‘1.4096-02 

12 

1.3226-02 

13  — 

1.4206-02 . 

14 

1.822E-02 

15 

1.5826-02 

16 

9.4346-03 

17 

8.5076-03 . 

18 

1.0216-02 

19 

1.1606-02 

20 

1.275E— 02 

~  21 

1.0946-02 

22 

9.5136-03 

23  . 

9.2296-03 

24 

8.0476-03 

25 

8.4026-03 

~ 3.4426-01 
_5. 3086-02 
2*6936-0 2 
2*1336-02 
'  1.2596-02 
7.6346-03 
1*3916-02 
1* 8416-02 
“1.3146-02‘ 
9.3436-03 
1. 0606-02 
1.2516-02 
1.3716-02 
1.4076-02 
‘1.2026-02" 
7.6096-03 
5.8846-03  ' 
8.0426-03 
1.0766-02  '* 
9.5736-03 
6.7266-03  ~ 
7.9786-03 
8.0736-03 
5.4876-03 
3.9426-03 


26 

1.1226-02 

4.2946-03 

27 

1.5786-02 

7.5136-0  3 

28 

1.4646-02 

8.8186-03 

29 

1. 0286-02  "*** 

6.1116-03 

30 

9.7406-03 

3.7486-03 

31 

9.4866-03 

‘  3.460E-0  3 

32 

7.3916-03 

4.6226-03 

33" 

4.7876-03 - 

6.8666-03 

34 

4.7156-03 

9.2876-03 

35 

5.8576-03  . 

“  8.5486-03 

36 

5.4626-03 

5.9226-03 

37  ‘ 

5.9526-03 

4.574E— 03 

38 

6.030E-03 

5.105E— 03 

39 

5.9296-03 

6.479E— 03  * 

40 

7.8576-03 

7.1226-03 

4l 

8.2366-03 

9.1776-03 

42 

8.3766-03 

9.0946-03 

43 

8.8186-03  . 

5.224E-03 

44 

8.3876-03 

4.352E-C3 

45 

8.8886-03 

5.8006-03 

46 

7.7106-03 

7.0966-03 

47 

6.6956-03 

6.508E-03 

48 

7.9196-03 

4.4336-03 

49 

7.5766-03  - 

5.1316-03 

50 

6.576E-03 

6.3726-03 

2.5496-01 

_ 1.4056-03 

1.7216-03 
5.117E-03 
2.1126-03 
1.121E— 03 
2.2856-03 
_ 1.4766-03 
-7,7046-05" 
-2.5526-03 
‘  -3.321E-03 
1.7036-04 
“  1.2116-03 

_ -2.8616—03 

-4.3946-03  ~ 
-3.0996-03 
”  -1. 7666-03 
1.5656-03 
2.8596-03 
1.7196-03 
2.1476-03  ~ 
3.0286-03 
1.7776-03 
-1.3316-03 
“-1. 3266-03 
_-4. 9096-04 
-3. 3916-03  “ 
-3.8456-03 
"  -1.158E-03 
-1.263E-03 
-1.4106-03 
-7.361E-04 
-5.0666-04 
-1.664E-04 
-9.511E-05 
1.991E-04 
-4. 9606-04  “ 
-1.631E— 03 
-1.1836-03 
-4.772E-04 
-1.598E-C3 
-3. 867E— 03 
-2.9586-03 
6.C95E-04 
1. 122E-03  - 
4.6866-04 
1.0396-03 
-3.506E-04 
-8.CC2E-C4 
3.Q21E— 04 


4.6476-02 

_  1.5326-02 

5.79SE— 03 
3.2696-03 
2.6956-03 
1.1556-03 
-2.079E-03 
_-3.354E-03 
-1.5316-03 
1.1346-03 
1.5916-03 
-4.963E-04 
— 2.805E-03 
-3.1056-03 
-1.178E-03 
3.777E-04 
.  -1.7606-04 
- -2.040E— 03 
“-3. 591 6-03 

_ —3. 181E— 03 

-5.3446-04 
9.8386-04 
1.7866-04 
-2.8616-05 
4.780E— 04 
1.0866-03 
7.7696-04 
4.450E-04 
7. 118E-04 
9.0566-04 
9.6776-04 
2.2716-04 
"-9. 2576-04 
.  -1.2366-03 
2.524E-04 
7.6856-04 
-3.2556-C4 
— 8.664E— 04 
-1.184E-03 
-3.805E-04 
1.1986-03 
1.2966-03 
8.8526-04 
2.8636-04 
-1.902E-C3 
-3.1816-03 
-1.7096-03 
-8.2326-04 
-8.7016-04 
0. 


■^•TV’'-’-’  %3sk,  ffr/.j|t;#(iB.)_i 


SPECTRAL  ESTIMATES  FOR 

ITEM  A  WITH  ITEM  B 


SILL  4  LINE  1C  2  TRACK  -40- 


p 

WXCP1 

0 

1 

7.6386-01 

2 

6.560E-02 

3 

1. 8876-02 

4 

2.1376-02 

‘  5 

2.4926-02 

6 

2.2796-02 

7 

2.0546-02 

8 

1.8736-02 

1.3146-02 

16 

9.4016-03 

iir 

7.4756-03 

12 

5.4276-03 

13 

6.8176-03 

14 

7.6196-03 

15 

6.5646-03 

16 

7.655E-C3 

17 

8.933E-03 

18 

9.7976-03 

19 

1.0096-02 

20 

7.1356-03 

21 

5.0746-03 

22 

7. 7396-03 

23 

1.0546-02 

24 

1.0216-02 

25 

7.9016-03 

26 

6.1086-03 

27 

7.2966-03 

28 

7.9356-03 

29 

5.5406-03 

30 

5.6736-03 

31 

8.5896-03 

32 

9.0756-03 

33 

8.5806-03 

34 

9.2476-03 

35 

8.5146-03 

36 

7. 5116-03 

37 

8.6006-03 

38 

8.6586-03 

‘  39 

_6. 569E-C3 

40 

5.433E-03 

41 

5.4286-03 

42 

5.5226-03 

43 

5.8966-03 

44 

6.5166-03 

45 

•6.1906-03 

46 

4.1466-03 

47 

3.1576-03 

48 

4.5016-03 

49 

5.707E-03 

50 

5.9366-03 

WY(P) 


1.158E  00 
1.0316-01 
4.0126-02 
3.356E-02 
1.6386-02 
8.262E-03 
9.3006—03 
1.1946-02 
1.4976-02 
1.5096-0  2 
1.5596-02 
1.4626-02 
1.1276-02 
1.1406-0  2 
1.0376-02 
7.4136-03 
8.0676-03 
8.9066-03 
7.7516-03 
6.901E— C3 
6.3286—0  3 
6.2116-03 
6.790E— 0  3 
6.711E-03 
5.9566-03 
4.5186-03 
5.857E— 03 
8. 1486-03 
5.7816-03 
4.2276-03 
'  5.588E-03 
5.1576-03 
3.878E-03 
5. 3076-03 
'  7.1796-03 
8.5786-03 
8.9056-03 
6.3846-03 
6.2816-03 
8.3346-03 
7.162E-03 
5.0306-03 
3.6656-03 
3.3096-03 
5.3646-03 
8.4886-03 
9.3926-03 
6.1926-03 
4.8376-03 
6.1346-03 


WC(P) 


8.0636-01 
5  •  867E-02 
8.7326-03 
5.033E-03 
3.1256-03 
1.0366-03 
1.786E— 03 
2.5746-03 
'  6.3656-04 
-1.3396-03 
-2.C216-Q3 
-1.4636-03 
2.8346-05 
2.2686-04 
~  8.6316-04 
1.5456-03 
-7.7996-04 
-1.3876-03 
-5.7876-05 
-1.2566-03 
— 2.3C46-03 
-8.4746-04 
7.669E-C4 
8. 139E-04 
1 .2906—03 
1.7906— C3 
1.443E-03 
5.361E— 04 
-8.6946-04 
-6.7026-04 
4. 455E-04 
1.6636-03 
'  2 . 8C3E— 03 

3.C426-03 
2.7686-03 
2.678E-03 
2.0526-03 
6.8446-04 
3.0206-04 
6.1536-04 
4.4446-04 
-9.6906-05 
1.7726-05 
4.4556-04 
-5.5006-04 
-1.8006-03 
-1.9256-03 
-1. 0666-03 
6.6946-05 
5.0736-04 


WQ{  P) 


2.4846-02 
— 1.207E— 02 
-7.6886-03 
2.771E-03 
2.9316-03 
-1.511E-03 
-2.4976-03 
1.405E-03 
4.6906-03 
2.402E-03 
— 2.849E— 04 
-1.4506-03 
-2.9466-03 
-1.3366-03 
8.3336-04 
3.2876-04 
1.0896-03 
3.145E-03 
3.5206-03 
2.503E— 03 
6.2296-04 
-1.5786-03 
— 2.783E— 03 
— 3.687E— 03 
-3.2946-03 
-1.5326-03 
-1.0726-03 
-3.8676-04 
8.9136-04 
1.144E— 03 
1. 105E-C4 
-4.1256-04 
-2.7866-0,4 
-1.332E-C3 
-2.4306-03 
-2.5736-03 
-i. 5766-03 
-5.8206-04 
9.971E-04 
2.6316-03 
1.716E— 03 
3.7866-04 
-4.4486-04 
-7.4526-04 
4.6096-04 
9.6446-04 
-3.8906-04 
-1.183E-03 
-5.0946-04 
0. 


SPECTRAL  ESTIMATES  FOR  SILL  1A  LINE  L  2  TRACK  -41 


ITEM  A  WITH 

ITEM  B 

p 

WX(P> 

WY(P) 

NC(P) 

MQ(P) 

u 

* 

1 

6.7356  00 

8.0706  00 

6.6076  00 

-2.6436  00 

2 

7*2106-01 

1.204E  OC 

7.1346-01 

-1.8586-01 

3 

1.890E-C1 

2.674E-01 

1.354E— 01 

-1.7326-02 

4 

9.0476-02 

6.703E-02 

4.1406-02 

. -1.6416-02 

5 

5.561E-02  ' 

2.6626-02 

2.0946-02 

-3.9086-03 

6 

4.214E-02  . 

1 .9296-02 

1.564E-02 

9.358E-C4 

7 

3.7466-02 

1.8886-02 

9.8386-03 

5.5696-04 

8 

2.9566-02 

2.5236-02 

5.5966-03 

-7.5496-05 

9 

2« 1936-02  ‘ 

3.507E-02 

"  9.0656-03 

-4.137E-03 

10 

2.631E-02 

3.1986-02 

9.4006-03 

-2.153E-03 

11 

2.6636-02  ' 

2.4386-0  2 

'*■  2.0136-03 

1.9696-03 

12 

1.7496-02 

2.873E-02 

-3.9686-03 

3.858E-03 

13 

. 1.8556—02 

"""  2.932E— 0 2 

-2.9446-03 

5.2676-03 

14 

2.5486-02 

2. 2446-02 

2.5176-03 

7.433E-03 

15 

2.7756-02 

2.8436-02 

6.5666-03 

8.5456-03 

16 

2.4936-02 

3.4756-02 

6.5736-03 

1.2786-03 

17 

1.8316-02 

2.8556-02 

2.6276-03 

-2.410E— 03 

18 

1.5096-02 

2.3916-02 

-3.2186-03 

1.145E— 03 

19 

1.4386-02 

'  2.1356-02 

-3.8106-03 

-8.0756-04 

20 

1.4496-02 

'1.9646—0  2 

-9.2096-04 

— 2.679E— 04 

21 

2.0856-02 

1.9506-02 

-4.3356-03 

2.817E— 03 

22 

2.775E-C2 

1.6716-02 

-7.5386-03 

-2.1006-03 

23 

2.6046-02 

"""*  1.4586-02 

— 4.962E— 03 

-4.0586-03 

24 

2.2836-02 

1.5896-02 

-2.1466-03 

-1.5196-03 

25 

~  1.9026-02  ‘ 

1.5326-02 

-4.3456-04 

-4.5566-03 

26 

1.1596-02 

1.4006-02 

-4.4356-04 

-3.757E-C3 

27 

1.0776-02 

1.7366-02 

5.9276-04 

2.51CE-03 

28 

1.651E— 02 

1.5826-02 

2.4796-03 

4.8186-03 

29 

1.8796-02 

1.0106^02 

2.4096-03 

3.4236-03 

30 

1.5026-02 

1.0866-02 

1.2296-03 

2. 3146-04 

31 

1.531E-C2 

1.1336-02 

-2.7536-03 

. -1.9526-03 

32 

2.181E-02 

1.6296-02 

-2.0716-03 

4. 557E-C3 

33 

1.914E-02 

2.8466-02 

4.641E-03 

7.9146-03 

34 

8.9866—03 

2. 826E-02 

5.5916-03 

2.5336-03 

35 

7.5536-03 

2.4196-02 

2.9766-03 

2. 7676-04 

36 

1.1036-02 

2.2976-02 

-1.1426-04 

1.2526-04 

37 

1.3686-02 

1.5456-02 

-2.7916-04 

2.8976-04 

38 

1.7296—02 

1.0946-02 

3.2416-03 

1.427E-03 

“  39  ‘ 

1.5976-02 

1.0996— 0  2 

2.937E— 03 

1.9366-03 

40 

9.981E-03 

8.0306-03 

1.2926-03 

1. 8666-03 

41 

7.972E-03 

7.5846-02 

1.7446-03 

2. 9656-03 

42 

1. 121E-C2 

1.2546-02 

9.4266-04 

4.684E-03 

43 

1.4146-02 

1.5486-02 

-1.5356-04 

3.7756-03 

44 

1.2026-02 

1.9166-02 

3.3736-03 

6.3076-04 

45  * 

1.0586-02 

2.4706-02 

“  7.8366-03 

-1.1446-03 

46 

1. 377E-C2 

2.0606-02 

4.6006-03 

5.4016-04 

47 

1.5286-02 

1.3086-02 

5.6176-04 

1.6376-03 

48 

1.4146-02 

1.7606-02 

1.5546-03 

-6.9916-04 

49 

1. 391E-02 

2.8386-02 

4.4C26-04 

-1.9506-03 

50 

1.432E-02 

3.3036-02 

-1.2896-03 

0. 

SILL  2A  LINE  1  TRK 


SPECTRAL  ESTIMATES  FOR 


ITEM  A 


p 

WX(P) 

0 

— 

1 

2.546E-01 

2  . 

2.043E-02 

3 

1.063E-02 

4 

7.964E-03 

5 

5.857E-03 

6 

4.677E-03 

7 

3.856E-03 

3 

3.767E-03 

9 

3.413E-03 

10 

2.494E-03 

11 

1.987E-03 

12 

1.957E-03 

13 

2.420E-03 

14 

2. 5536-03 

15 

1.692E-03 

16 

1. 325E-03 

17 

1.815E— 03 

18 

1.821E-03 

19 

1.444E-03 

20 

1.327E-03 

21 

1.299E— 03 

22 

1.251E— 03 

23 

1.000E-03 

24 

8.265E-04 

25 

1.087E-03 

26 

1.105E-03 

27 

9.580E-04 

28 

8.899E-04 

29 

6.888E-04 

30 

6.761E-04 

31 

7.910E-04 

32 

7.827E-04 

33 

8.633E-04 

34 

9.906E-04 

35 

1.120E-03 

36 

1.045E-03 

37 

7.976E-04 

38 

8.580E— 04 

39 

7.903E-04 

40 

4.961E-04 

41 

4.901E-04 

42 

6.271E-04 

43 

6.490E-04 

44 

6.411E-04 

45 

5.758E-04 

46 

3.677E-04 

47 

3. 079E-04 

48 

3. 832E-04 

49 

5.323E-04 

50 

7.716E-04 

SPECTRAL  ESTJ 
ITEM  A 


0 

1 

2 

3 

4 

5 

6 
7 
3 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 


MATES  FOR 


HOOD  1  LINE  N 


WX(P) 


3.876E  00 
1.138E  00 
3«4)70E— 01 
8.701E-02 


2.3| 

l.ll 


71E-02 

04E-02 


1.049E— 02 

1.139E-02 

1.304E-02 

1.558E-02 

1.2^5E-02 

7.4|19£— 03 

5.910E-03 

6.767E-03 

6.680E-03 

5.441E-03 

4«4[73E-03 

5.409E-03 

6.756E-03 

6.3 75E-03 

4.5£46-03 

3.4B2E-03 

3.8B0E-03 


25 

4.8 

|ri£-03 

26 

6.3 

77E-03 

27 

6.2 

95E-03 

28 

4.9 

L7E-03 

29 

4.2^96-03 

30 

3.631E-03 

31 

3.204E-03 

32 

4.1^76-03 

33 

5.626E-03 

34 

7.015E-03 

35 

8.010E-03 

36 

5.707E-03 

37 

4.472E-03 

33 

5.973E-03 

39 

5.967E— 03 

40 

7.327E-03 

41 

7.  71 

6E-03 

42 

5.  36 

56-03 

43 

4. 556E-03 

44 

5.015E-03 

45 

5.803E-03 

46 

5.330E-03 

47 

4.188E-03 

48 

3.228E-03 

49 

3.44 

5E-03 

50 

3.55 

AF-03 

-43 


^  --.  ■■■*..*  , ^  _ _ _  _ 

SPECTRAL  ESTIMATES  FUR 

ITEM  A 

P  W X(P) 


HOOD  2  LINE  0  -44. 


0 


1 

5.103E-01 

2 

6.514E-02 

3 

2.797E-C2 

4 

1.971E-02 

5 

1.201E-02 

6 

6.580E-03 

7 

6.232E-03 

8 

7.6S4E-03 

9 

8.986E-03 

10 

1.129E-02 

11 

1.423E-02 

12 

1./63E-02 

13 

1.920E-02 

14 

1.593E-02 

15 

1.072E-02 

16 

6.966E-03 

17 

4.266E-03 

IS 

2.297E-03 

19 

1.824E— 03 

20 

2.131E-03 

21 

2.761E-03 

22 

3.380E-03 

23 

3.214E-03 

24 

2.596E-03 

25 

1.8?'  -03 

26 

l«70bt  03 

27 

2.057E-03 

28 

1.663E-03 

29 

1.251E— 03 

30 

1.117E-03 

31 

8.391E— 04 

32 

8.855E-04 

33 

1.350E-03 

34 

1.879E-03 

35 

2.015E-03 

36 

2.085E-03 

37 

2.341E-03 

38 

2.101E-03 

39 

1.915E-03 

40 

2.167E-03 

41 

2. 101E-03 

42 

1. 715E-03 

43 

1.667E-03 

44 

1.838E-03 

45 

1.629E-03 

46 

1.423E-03 

47 

1.314E-03 

48 

1.092E-03 

49 

1.206E-03 

50 

1.415E-03 

\ 

! 


SPECTRAL  ESTIMATES  FOR  POLK  1  LINE  P  -45- 


ITEM  A 


p 

WXIP) 

0 

1 

5.966E  00 

2 

7.1116—01 

3 

4.693E-01 

4 

2.892E-01 

5 

1.555E-01 

6 

9. 1706-02 

7 

6.343E-02 

8 

5.137E-02 

9 

5.256E-02 

10 

5.533E-02 

11 

5.101E-02 

12 

5.379E-02 

13 

5.920E-02 

14 

5.263E-02 

15 

4.841E-02 

16 

5.835E-02 

17 

6.054E-02 

18 

4.147E-02 

T9 

3. 2816— 02 

20 

4.090E-02 

21 

4.1876-02 

22 

3.468E-02 

23 

2.616E— 02 

24 

2.2286—02 

25 

1.787E-02 

26 

1. 110E-02 

27 

9.4686-03 

28 

1.080E-02 

,29 

1.1916-02 

30 

1.3656-02 

31 

1.614E-02 

32 

1.815E-02 

33 

1.755E-02 

34 

1.668E-02 

35 

1.458E-02 

36 

9.7636-03 

37 

6.975E-03 

38 

7.8006-03 

39 

8.930E-03 

40 

7.545E-03 

41 

5.7106-03 

42 

5.646E-03 

43 

5.5836-03 

44 

5.628E-03 

45 

5.972E-03 

46 

4.9536-03 

47 

3.2656-03 

48 

3.684E-03 

49 

6.2666—03 

50 

7. 77 IE- 03 

ni 


SPECTRAL  ESTIMATES  FOR  POLK  2  LINE  Q 

item  a 


p 

WX(P> 

0 

.  ---- 

1 

8.006E  01 

_  2 

9.443E  00 

3 

3.610E  00 

4 

1.473E  00 

5 

1.169E  00 

6 

1.202E  00 

7 

8.202E-01 

8 

4.226E—0I 

9 

2.274E-01 

10 

1.171E-01 

11 

6.953E— 02 

.  12 

8.722E-02 

13 

1.097E-01 

14 

8. 112E-02 

15 

8.318E-02 

16 

1.570E-01 

17 

1.518E-01 

18 

7.988E-0 2 

19 

9.450E-02 

20 

1.282Er01 

21 

8.606E-02 

22 

5.780E— 02 

23 

5.860E-02 

24 

4.296E-02 

25 

4.932E-02 

26 

6.210E— 02 

27 

4.751E— 02 

28 

3.146E-02 

29 

3.371E-02 

30 

4.429E-02 

31 

5.382E-02 

32 

4.898E-02 

33 

3.489E-02 

34 

2.461E-02 

35 

1.987E-02 

36 

2.003E-02 

37 

1.984E-02 

38 

1.959E-02 

39 

1.573E-02 

40 

8.489E-03 

41 

6.527E-03 

42 

1.219E-02 

43 

2.114E-02 

44 

2. 036E-02 

45 

1. 108E-02 

46 

8.019E-03 

47 

1.093E-02 

48 

1.153E-02 

49 

1. 114E-02 

50 

1.208E-02 

SPECTRAL  ESTIMATES  FOR 


ITEM  A 


HXCFJ 


BENNI NS  J^LINE  R 


'  * '  7  : '  i  v  *- 


1 

_ 2 _ 

3 

4. 

5 

6 
7 

_  8 _ 

9 

10 

11  :  " 
12 

13 

14 
~  15 

16 

17  "" 

18 

19 

20  _ 

“21 

22 

23 

24  : 

‘25'"' 
_26 _ 

27 

28 

29 

30 

31 
-32 

33 

.34..... 

35 

.  36  _  _ 

37 

38  _ 

39  ' 

40 

41  " 

42  .. 

43 

44  _ 

45 

46 

47 
40 

49 

50 


6.1356 
1.328E 
6.7366- 
2.855E- 
1.328E— 
1.204E- 
9.200E- 
4.219E- 
3.108E- 
4.6396- 
4.549E- 
3.441E- 
'  3.1596- 
3.6006— 
3.6296- 
2.554E- 
1.4896- 
1.1246- 
1.2566- 
1.6476- 
1.791E- 
1.7186- 
1.6266- 
1.3596- 
1.0706- 
8.959E- 
8.736E- 
1.1796- 
1.3326- 
1.185E- 
1.0206- 
7.746E- 
9.2976- 
1.425E- 
1.4586- 
1.241E- 
9.9526- 
7.3306- 
5.1366- 
3.9676- 
5.5486- 
8.0756- 
1.010E— 
1.301E- 
1.433E- 
1.1696- 
9.4616- 
7.7176- 
5.6036- 
4.974E- 


f 

I 


„_$.PECTRAl  _£ST  I  MATES .  FOR 

_  :  ITEM  A  _ ' 

F  NX  IP)  J _ 

„  ®— -*■  u  l45g  00  ^ 

_ ; _ 2j _ 1.566E-01 _ _ 

3  1.145E-01 

_ 4  _ 1  •  44 9 E-  0 1 _ 

5  1.522E-01 

6  "  1.102E-01  _  . 

7  6.686E-02  , 

_  _  4.708E-02  ^ _ 

9  .  3.425E-02 

10  2.275E-02  • 


Hi/ 

1.723E-02 

12 

1.556E-02 

13 

1.742E-02  . 

14  . 

1.936E-02 

15 

1.369E-02  , 

16  .  . 

8.478E-03  ^ 

17 

1.186E-02  , 

18  / 

1.583E-02 

19 

1.373E-02  , 

20 

1.160E-02  __ 

21  / 

1.476E-02  . 

22 

1.692E-02  i 

23  , 

1.650E-02 

24 

1.659E-02 

25 

1.446E-02  , 

26 

?,.048E-02 

27 

1.039E-02 

28  „r„. 

1.232E-02  .. 

29 

1.083E-02 

30 

8.285E-03  „ 

31 

8.183E-03 

32 

7.564E-03 

33 

9.945E-03  , 

34 

1.503E-02 

35 

1.285E-02  , 

36  . 

7.019E-03  . 

37 

4.865E-03  , 

39 

•5.648E-03 

39  .  ... 

8.056E-03 

9.917E-03  „ 

41 

8.673S-03  „ 

42 

6«083c-03  j . 

43 

5.560E-03 

44 

5.904E-03 

45 

5.437E-03 

46 

5.589E-03  . 

47 

6.196E-03 

40 . 

'6*  36  20-03  . 

49  5.336C-03 


BENNIN6  2. LINE  $ 


s.  Jr  'nUJ"V '*'«<•- 


i 


SPECTRAL  ESTIMATES  FOR  MCCLELLAN  1  LINE  T  -48- 


ITEM  A 


p 

WXCP) 

c 

— — 

1 

1.711E  00 

2 

5.754E-01  ... 

3 

2*937E— 01 

4 

1.842E-01 

5 

1.276E-01 

6 

7*948E— 02 

7 

5«577E— 02 

8 

4.659E-02 

9 

4.523E-02 

10 

3.830E-02 

11 

2.450E-02 

12 

2.144E-02 

13 

2.398E-02 

14 

1.843E-02 

15 

1.024E-02 

16 

8.966E-03 

17 

1.365E-02 

18 

1.603E-02 

19 

1.265E-02 

2C 

9.566E-03 

21 

9.606E-03 

22 

!• 162E—02 

23 

1.255E-02 

24 

1.291E-02 

25 

1« 169E— 02 

26 

7.310E-03 

27 

7.545E-03 

28 

1.297E-02 

29 

1.747E-02 

30 

1.617E-02 

31 

1. 212E-02 

32 

8.836E-03 

33 

6«726E— 03 

34 

6 *24 IE— 03 

35 

7.803E-03 

36 

1.104E-02 

37 

1.232E-02 

38 

1.109E-02 

39 

9.863E-03 

40 

8.742E-03 

41 

1.050E-02 

42 

1.259E-02 

43 

1.031E-02 

44 

9.272E-03 

45 

9.000E-03 

46 

7.964E-03 

47 

8.902E-03 

48 

9.391E-03 

49 

8.535E-03 

5C 

7*907£-03 

SPECTRAL  ESTIMATES  FOR 


MCCLELLAN  2  LINE 
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AREA  SPECTRA 


Again  our  methods  of  obtaining  p.s.d.  estimates  of 
area  data  are  described  In  [1].  The  computer  program 
used  will  be  described  In  the  section  on  computations. 

In  this  section,  we  present  the  spectral  estimates 
In  numerical  form  and  contour  graphs.  The  raw  data  are 
too  voluminous  to  present  in  this  report;  persons  wish¬ 
ing  to  use  this  data  may  obtain  it  at  cost  In  the  more 
convenient  form  of  IBM  cards  or  magnetic  tapes  from 
MASC. 

Because  of  symmetry,  only  half  of  the  spectral  esti¬ 
mates  are  presented  In  tabular  form.  Somewhat  more  com¬ 
plete  contour  graphs  are  given.  The  contour  heights  are 
In  powers  of  5  corresponding  to  the  seml-logarlthmlc 
plots  of  the  line  spectra.  On  each  contour  diagram  Is 
an  arrow  indicating  the  direction  North. 

The  statistical  accuracy  of  the  spectral  estimates 
may  be  presented  as  on  p  72  of  [1].  Each  estimate  has 
32  degrees  of  freedom.  A  95Z  confidence  Interval  may 
be  given  by  the  factors  .648,  1.740w  [We  are  95J  con¬ 
fident  that  a  true  spectral  value  Is  between  .648  times 
the  estimate  and  1.740  times  Its  estimate.] 
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0.0494407 

0.0460097 

0.0410304 

-3 

0.0659439 

0.0573596 

0.0488996 

i  -2 

0.0919362 

0.1062143 

0.0958538 

-1 

0.3155221 

0.3246462 

0.3354494 

0 

0.5477226 

0.4725607 

0.5111475 

I  1 

0.3075510 

0.280958 2 

0.2826260 

1  '  2 

0.0931707 

0.1142454 

0.0945172 

|  3 

0.0620384 

Q. 0602905 

0.0429307 

4 

0.0376946 

0.0298909 

0.0238513 

5 

0.0296627 

0.0265874 

0.0245691 

6 

0.0347422 

0.0304845 

0.0263233 

7 

0.0375168 

0.0315281 

0.0252664 

a 

0.0328973 

0.0305252 

0.0315638 

9 

0.0280882 

0.0276388 

0.0317478 

10 

0.0268868 

0.0253586 

0.0260205 

11 

0.0271181 

0.0243949 

0.0241802 

12 

0.0279903 

0.0279021 

0.0297137 

13 

0.0214373 

0.0219197 

0.0226373 

14 

0.0173777 

0.0193397 

0.0225023 

15 

0.0186926 

0.0228444 

0. 0245694 

16 

0.0150734 

0.0204429 

0.0196998 

17 

0.0161044 

0.0196422 

0.0235398 

18 

0.0188831 

0.0204996 

0.0280960 

19 

0.0172860 

0.0176335 

0.0207658 

20 

0.0159938 

0.0161634 

0.0188156 

•  •» 

I 
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COLUMNS  16  TO 

18 

16 

17 

18 

* 

-20 

0.0175210 

0.0141239 

0.0150464 

-19 

0.0159988 

0.0144360 

0.0166053 

* 

-18 

0.0167070 

0.0159709 

0.0194154 

-17 

0.0196668 

“  0.0165106' 

0.0178439 

-16 

0.0209756 

0.0172503 

0.0186742 

* 

-15 

0.0239890 

0.0214877 

0.0217651 

j 

-14 

0.0247584 

0.0208252 

0.0186370 

j 

-13 

0.0223732 

0.0213959 

0.0203281 

! 

-12 

0.0209229 

0.0226293 

0.0261093 

1 

-11 

0.0243640 

0.0268335 

0.0296806 

i 

-10 

0.0294413 

0.0278359 

0.0260897 

-9 

0.0274384  ~ 

r0.027?405 

"  0.0282074  ' 

; 

-8 

0.0270444 

0.0259220 

0.0289039 

-7 

0.0340487 

0.0371019 

0.0429154 

■a 

-6 

0.0383115 

0.0449601 

0.0498114 

? 

-5 

0.0406082 

0.0469148 

0.0509506 

j 

-4 

0.0407289 

0.0396328 

0.0433115 

! 

-3 

0.0448129 

0.0395367 

0.0378795 

i 

4 

& 

_  . _  _ _  _  . .  i 

-2 

0.0831960 

0.0868085 

0.0733807 

-1 

0.3073666 

0.3452572 

0.3738297 

* 

0 

0.4633564 

0. 4906546 

0.6155693 

w 

I 

1 

0.2303791 

0.2472107 

0.3557941 

2 

0.0651278 

0.0794220 

0.1115352 

| 

3 

0.0373373 

0.0405093 

0.0441264 

4 

0.0277572 

0.0315612 

0.0339950 

r  ‘  5 

0.0264499 

0.0265586 

0.0255295 

4 

~  6 

0.0279019 

0.0299479 

0.0304182 

7 

0.0289191 

0.0302710 

0.0325101 

1 

8 

0.0357182 

0.0367567 

0.0342234 

i 

i 

;  9 

0.0304774 

0.0286684 

0.0248205 

2 

j 

-  10 

0.0292091 

0.0285102 

0.0188145 

] 

11 

0.0305088 

0.0298469 

0.0204559 

J 

12 

0.0324184 

0.0290934 

0.0243034 

13 

0.02646C  0 

0.0233542 

0.0232352 

■ 

14 

0.0254737 

0.0238434 

0.0232587 

| 

15 

0.0242747 

0.0211017 

0.0203917 

! 

16 

0.0161775 

0.0169190 

0.0197144 

» 

j 

17 

0.0214768 

0.0192246 

0.0200203 

! 

18 

0.0252437 

0.0191662 

0.0205691 

*  ! 

19 

0.0182629 

0.Q162892 

0.0174786 

20 

0.0167277 

0.0180698 

0.0191875 

RILEY 


RILEY 


GRUBER 


f  cpf 


CARSON 


-.20 


SILL 


f  cpf 


McClellan 


McClellan 


ABERDEEN 


ABERDEEN 


YUMA 


YUMA  2 


COMPUTATIONS 


The  computations  for  this  report  are  done  on  two 
programs  for  the  IBM  709*1  computer.  These  programs  will 
be  described  in  considerable  detail  in  Appendices  B  and 
C. 

In  Appendix  B#  is  the  co-spectral  program  for  process 
lng  data  from  one  track  or  two  parallel  tracks.  Tnls 
program  was  modified  from  an  existing  program  by  adding 
to  it  a  routine  for  smoothing  the  input  data  and  correct¬ 
ing  the  resulting  spectrum.  Other  changes  of  a  minor 
nature  were  also  made. 

The  program  described  in  Appendix  C  was  made  espec¬ 
ially  for  processing  area  data;  it  is  new. 

Both  programs  are  available  on  IBM  cards  at  Midwest 
Applied  Science  Corp.  at  cost. 


-216- 


CONCLUSIONS 


There  were  no  surprises  in  the  new  line  data.  Most 
of  the  features  presented  in  the  nri.d.'s  of  our  previous 
report  [1]  are  again  to  be  obs^rtfed  in  the  new  p.s.d.'s. 
However,  the  bumps  near  the  2?f  ft.  wave  lengths  noticeable 
In  the  Knox  and  Aberdeen  p.s.d.’s  are  not  observed  again. 

We  conjectured  in  [1]  that  these  bumps  were  due  to  periodic 
components  observable  in  the  ground  profile.  These  could 
have  been  caused  by  repeated  use  of  the  ground  by  vehicles. 

The  spectra  trail  off  to  levels  associated  with  meas¬ 
urement  error,  as  noted  in  [1].  Thus,  power  in  the  high  fre¬ 
quency  components  was  generally  low.  This  suggests  that  a 
longer  measurement  interval  could  have  been  used.  But  we  are 
constrained  to  note  that  in  one  vehicle  test  at  Port  Sill 
high  frequency  vibration  present  in  the  steering  system  tend¬ 
ed  to  limit  speed,  indicating  that  in  certain  situations  clos 
er  spacing  of  data  points  may  be  required;  this  point  will  be 
commented  upon  again  in  our  report  on  vehicle  vibration  tests 

Since  Interpretation  Is  reserved  for  the  second  report, 
we  shall  now  confine  attention  to  conclusions  concerning  data 
acquisition  and  processing. 

Surveying  methods  present  a  very  simple  technique  for 
acquiring  elevation  data.  Equipment  is  easy  to  rent  and  use. 
It  is  easy  to  train  unskilled  personnel  in  their  use  for  this 
purpose.  The  data  from  Las  Vegas  were  taken  in  one  day  by 
persons  who  had  not  previously  made  this  type  of  measurement. 

One  Is  tempted  to  draw  general  conclusions  from  a 
visual  inspection  of  the  data  presented.  If,  however, 
such  remarks  are  to  be  meaningful,  they  must  be  made  keep¬ 
ing  in  mind  the  purpose  of  the  study.  We  cannot,  for  ex¬ 
ample,  assert  that  the  accuracy  of  our  estimates  is  good 
until  we  know  what  accuracy  is  necessary.  We  cannot  make 
up  curves  to  fit  the  results  unless  we  know  what  part  of 
the  results  must  require  the. closest  fit.  The-  Internal  con¬ 
sistency  of  the  results — the  comparison  of  the  line  data 
with  the  area  data  at  each  site,  for  example — requires  some 
analytic  work  and  perhaps  more  computation.  We  will  delay 
such  conclusions  to  the  second  report.  • 


The  spectral  estimation  procedure  was  designed  with 
conventional  vehicle  sizes  and  speeds  In  mind.  As  men¬ 
tioned  elsewhere.  It  Is  conceivable  that  reprocessing 
of  the  data  may  be  required  for  unconventional  vehicle 
types.  This  point  will  be  studied  in  due  course* 
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RECOMMENDATIONS 


1.  Pop  measurement  bf  ground  roughness,  we  recommend  that 
the  surveys  be  conducted  with  self-leveling  levels  and 
self-zeroing  rods.  Experienced  personnel  are  not  re¬ 
quired. 

2.  The  format  on  which  the  survey  data  is  entered  into 
books  is,  on  thei surface,  a  small  matter.  However, 
in  recording  it  bn  IBM  cards  for  program  input,  it 
may  be  a  matter  of  considerable  expense  if  the  orig¬ 
inal  data  are  not  in  a  suitable  format.  Recommenda¬ 
tions  for  recording  data  are  presented  with  the  pro¬ 
gram  descriptions  in  this  report. 

3.  A  pre-program  must  be  used  to  detect  outlying  data 
points.  The  anomalous  points  must  be  compared  with 
original  books  and  changed  1*  there  is  sufficient 
reason.  Unless  errors  of  th'  type  are  eliminated, 
p.s.d.  estimates)  will  have  ptcullar  features  which  are 
not  easily  explained.  Consultation  with  survey  per¬ 
sonnel  is  usually  desirable  to  check  whether  outlying 
values  are  errors. 

j 

4.  An  extended  program  to  obtain  a  larger  atlas  of  ground 
roughness  measurements  does  not  seem  justified  on  the 
basis  of  the  similarity  between  the  results  already  ob¬ 
tained,  Specific  experiments  or  vehicle  trials  may, 
however,  require  surveys. 


We  recor.jnend  this  atlas  be  used  by  those  interested 
in  aspects  of  ground  roughness  which  can  be  measured 
by  p.s.d.' s.  These  ; .s.d.'s  are  from  many  different 
areas.  They  na^  be  combined  into  an  average  for  some 


purposes  or  their  uifferences  may  be  noted.  In  any 
case,  a  quantitative  use^trf~these  results  awaits 
further  development.  (  \ 


REVIEWED: 


RONALD  A.  LISTON 
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APPENDIX  A 


List  of  Equipment  Used  by  Surveying  Crew 


dozen 

lb. 

rolls 

rolls 

dozen 


sets 


skeins 

sets 

set 


Blue  Crayons 
Stake  Tacks 
Red  Plastic  Tape 
Orange  Plastic  Tape  . 

Pencils 

Surveyors  Bags 

100  ft.  Engineer  Tapes  (Post  Cat.  No.  3141H) 
100  ft.  Lufkin  3/8"  Steel  Tapes 

(Post  Cat.  No.  3245K) 

100  ft.  Cloth  Tapes  (Post  Cat.  No.  3200K) 

Range  Poles 
14  In.  Marking  Pins 
Marking  Pin  Rings 
8  oz.  Plumb  Bobs 

Plumb  Bob  Sheaths  /' 

Plumb  Bob  Cord 
Level  Book  Leaves 
Cross  Section  Leaves 
Binders 

Prismatic  Compass 

Lenker  Elevation  Leveling  Rods 

Repeating  Theodolites  K&E  Cat.  No.  730050 

Zeiss  Self  Leveling  Levels  K&E  Cat.  No.  750020 

Auto  Top  Carrier  and  Cover 

Shovels 

Machettes 

Thermos  Jugs 

Folding  Stools 

Canvas  Bags 

Hatchet 

Snake  Bite  Kits 
Auto  First  Aid  Kit 
Camera 

Station  Wagon 

Surveying  Stakes  and  other  expendible 
items  purchased  in  field. 
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APPENDIX  B 

Linear  Power  Spectral  Penalty  Program 
'PSLINE* 


1.  Purpose:  'PSLINE'  la  designed  to  accept  one  or  more  sets 

of  data,  each  of  which  consists  of  one  to  five  parallel 
tracks.  The  tracks  are  represented  by  aa  many  as  one 
thousand,  equally  spaced  points.  Control  parameters 
have  been  provided  to  allow  flexibility  in  specifying 
filtering  values,  output  requirements  and  track  selec¬ 
tion. 

' PSLINE '  provides  the  basic  operation  of  computing 
the  covariance  functions,  raw  spectra  and  spectral 
estimates  for  one  track  or  a  pair  of  tracks. 


2.  Machine  Requirements ;  The  'PSLINE*  program  exists  as  a 
(IBM)  FORTRAN  iV source  program  consisting  of  several 
subroutines.  It  Is  Intended  to  be  run  on  the  IBM  7094 
computer  under  the  IBSYS  operating  system.  Deviation 
from  these  conditions  may  involve  slight  modification 
of  the  source  program.  (See  Supplement  1  for  detailed 
description  of  program  logic.) 


3.  Card  Preparation: 

CONTROL  CARD 

COLS  1-6 
7  -  8 


9 

10  -  12 
1.3 

14  -  15 


1LINEP 

Input  tape  number. 

¥5,  blank  -  all  Input  is  taken, 
from  tape  05. 

XX  -  TITLE  card,  track  data 
only  are  taken  from  tape  XX, 

'other  Information  from  tape  05* 

X  -  number  of  tracks. 

1  <  X  <  5 

XXX  -  number  of  lags. 

1 

XX  -  Number  of  filtering  coefficients. 
1  <  XX  <  20 


16  -  21  XXXXXX  -  distance  between  points 
within  a  track.  Decimal  point 
must  be  punched. 

22  0,  blank  -  do  not  print  raw  data. 
1  -  do  print  raw  data. 

23  0,  blank  -  do  not  print  smoothed 

data 

1  -  do  print  smoothed  data. 

24  >25  XX  -  number  of  selection  cards. 

26  -  80  Blank 


SELECTION  CARDS 
Type  1: 

COLS  1-7  ISELECT 

8  0  -  Print  covariance  functions 

1  -  Do  not  print. 

9  0  -  Print  raw  spectrum. 

1  -  Do  not  print. 

10  0  -  Print  spectral  estimates. 

1  -  Do  not  print. 

11  -  16  bbbbY  -  Track  No.  1 

17  -  22  bbbbb,  bbbbY  -  Track  No.  2 

In  both  cases,  if  Y  is  present 
1  <  Y  <  X,  where  X  i3  the  digit 
'  in~Col7  9  Of  the  CONTROL  card. 

If  Track  No.  2  is  present,  then 
the  pair  of  tracks  (No.  2  and 
No.  1)  will  be  processed. 

If  Track  No.  2  is  not  present, 
then  only  the  single  track 
specified  by  Track  No.  1  will 
be  processed. 

Type  2: 

COLS  1-7  $SELALP 

8,9,10  As  in  Type  1. 

All  possible  pairs  of  tracks 
are  processed. 
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Type  3i 

COLS  1-7  ISELALS 

8,9 ,10  As  in  Type  1* 

All  tracks  are  processed  Individually. 

KEY  CARD 

COLS  1-72  Blank 

73  -  80  An  eight  digit  number  which  appears 
on  the  TITLE  card  of  the  partic¬ 
ular  track  data  which  la  to  be 
processed. 

This  card  Is  present  to  allow  processing  of  arbi¬ 
trarily  positioned  data  sets  on  a  tape  other  than 
Tape  5.  However,  regardless  of  whether  the  track 
data  Is  on  Tape  5  or  on  another  tape,  the  KEY  card 
number  and  the  TITLE  card  number  for  the  track 
data  must  be  present.  If  more  than  one  set  of 
track  data  Is  to  be  processed,  the  sets  must  be  In 
the  same  order  as  the  KEY  cards  are  encountered, 
l.e.,  PSLINE  will  not  rewind  the  Input  tape  If 
the  correct  TITLE  card  has  not  encountered  when 
the  end  of  the  tape  Is  reached. 


TITLE  CARD 

COLS  1-6  Blank 

7-36  Any  text  which  identified  this 

track  data  set.  This  title  will 
appear  at  the  top  of  each  page 
of  output  for  this  track  data 
set . 

37  -  72  Blank 

73  -  80  Eight  digit  number  which  will 

identify  this  track  data  set 
given  the  number  on  the  KEY 
card. 


TRACK  DATA  A  set  of  track  data  consists  of 

1  or  more  tracks,  all  of  which 
consist  of  the  same  number  of 
points.  Let  T(I,J)  denote  the 
Ith  point  of  the  Jth  track.  Sup¬ 
pose  there  are  n  tracks ,  then 
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T(I,1) ,  T(I,2),  •••  T(I,n)  is 
called  the  Ith  cross-section 
for  this  set  of  track  data. 

Track  data  is  punched  into  cards 
by  cross-sections  and  is  read 
by  'PSLINS*  by  cross-sections. 

A  FORMAT  card  must  be  constructed 
to  indicate  how  the  cross-section 
values  have  been  punched  on 
-  cards.  In  addition  the  FORMAT 
card  must  allow  for  a  one  charac¬ 
ter,  alphameric  field  to  be  read 
after  the  last  number  of  each 
cross-section.  The  use  of  this 
--  character  will  be  discussed  later. 

Example  Track  14.0,4.2*4.3 
Track  2  3*1,  3.2,  3-3 


This  could  be  punched  as  follows: 


I 

V 

Card  1 

COLS 

1-3 

4.0 

jf 

i 

4-6 

3.1 

\ 

f 

Card  2 

COLS 

1-3 

4.2 

4 

4-6 

3.2 

Card  3 

COLS 

1-3 

4.3 

4-6 

3.3 

Card  4 

COLS 

1-3 

Blank 

4-6 

Blank 

7 

*  or  %  or  / 

And  the  associated  FORMAT  card  would  be 
(2F3.0,  Al) 


NOTE:  Columns  73-80  should  not  be  used  to  contain  track 
data  or  the  one  character  field. 


TERMINATION  CARDS 


In  order  to  eliminate  the  necessity  of  counting 
the  number  of  points/track  in  a  track  data  set, 
PSLINE  looks  for  a  special  cross-section  in  which 
the  one-character  field  is  not  blank  and  is  one 
of  *,  $  or  /.  When  a  cross-section  of  this  type 
Is  found  the  following  actions  are  taken: 

1.  The  numeric  values  read  are  NOT  Included 
in  the  true  track  data. 

2.  It  is  assumed  that  the  last  cross-section 
read  is  the  last  cross-section  for  this 
track  data  set. 

3.  The  particular  character  (*,  $,  or/) 
found  determines  how  the  NEXT  input  set 
is  to  be  processed: 

*  indicates  that  there  are  no  more 
input  sets. 

/  Indicates  that  the  next  input  set 
consists  of  only 

KEY  card 
TITLE  card 

Track  Data  (including  TERMIN¬ 
ATION  card) 

This  allows  one  to  process  many  track  data 
sets  with  the  parameters  found  on  one 
CONTROL  card.  Note  that  this  presupposes 
that  filtering  specifications,  number  of 
tracks ,  number  of  SELECTION  cards  and 
track  data  format  are  identical  for  the 
next  data  set. 

$  indicates  that  a  new  control  card  is 
to  be  read  with  the  next  input  set. 

I. el,  the  next  Input  set  will  consist 
of 

CONTROL  Card 

Filtering  Coefficient  Format  Card 
,  Filtering  Coefficients 

Track  Data  Format  Card 
KEY  Card 
TITLE  Card 


-  B-6  - 


Track  Data  (including  TERM¬ 
INATION  Card) 
SELECTION  Card(s) 


PORMAT  CARD 

The  variable  format  technique  is  employed 
by  PSLINE  to  allow  flexible  data  card  format. 

A  FORMAT  card  for  PSLINE  consists  of  a  standard 
FORTRAN  FORMAT  statement,  with  the  word  FORMAT 
deleted,  punched  free  form,  into  columns  1-72  of 
a  card. 

Filtering  coefficient  format  must  contain, 
only  F  or  E-type  conversions,  track  data  formats 
must  specify  one  more  field  than  the  number  of 
tracks  with  the  last  of  the  fields  Al. 

Consult  IBM  FORTRAN  IV  programming  (7090/9*0 
manuals  for  detailed  description  of  FORMAT  state¬ 
ments. 


4.  Deck  Preparation:  The  following  is  an  example  of  a  deck 

prepared  using  the  various  termination  character  options. 
Bracketed  cards  would  appear  on  another  tape  (other 
than  the  system  input  tape)  if  the  CONTROL  card  speci¬ 
fied  so. 

$J0B  card 

Installation  $ID  card  > 

$ EXECUTE  IB JOB 
IIBJOB 

Source  Program  Deck 
or 

Binary  Object  Deck 
$DATA  card 
CONTROL  card 

Filtering  Coefficient  FORMAT  Card 
Filtering  Coefficients 
Track  Data  FORMAT  Card 
KEY  Card 
TITLE  Card 

Track  Data  (TERMINATION  Card  with  /) 


SELECTION  Card( s ) 

KEY  Card 

Track  Data  (TERMINATION  Card  with  $) 
SELECTION  Card(s) 

CONTROL  Card 

Filtering  Coefficient  FORMAT  Card 
Filtering  Coefficient 
Track  Data  FORMAT  Card 
KEY  Card 
TITLE  Card 

Track  Data  (TERMINATION  Card  with  *) 
SELECTION  Card(s) 

END-OF-FILE  Card 
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SUPPLEMENT  1 

to  * 

LINEAR  POWER  SPECTRAL  DENSITY  PROGRAM 


A.  General  Program  Logic 

The  PSLINE  program  consists  of  a  main  program  (deck 
name  RVR...)  and  five  subroutines  (entry  points  -  INPUT, 
OUTPUT,  FILDEC,  COQUAD,  and  CORRET). 

1.  INPUT  -  The  INPUT  subroutine  reads  all  input 
cards  except  selection  cards  and  sets  parameters 
to  control  the  processing  of  the  track  data. 

Upon  reading  the  filtering  coefficients,  it  calls 
subroutine  CORRET  to  compute  the  actual  filtering 
values. 

2.  OUTPUT  -  The  OUTPUT  subroutine  is  called  by 
the  main  program,  RVR...,  to  print  out  raw  and/or 
smoothed  track  data,  one  track  per  call. 

3.  FILDEC  -  The  FILDEC  subroutine  is  called  by 
the  main  program,  RVR...,  to  perform  the  filter 
transformation  on  the  track  data. 

4.  CORRET  -  The  CORRET  subroutine  is  called  by 
the  INPUT  routine  to  compute  the  actual  filtering 
values  from  the  filtering  coefficients  read  from 
cards , 

5-.  COQUAD  -  The  COQUAD  subroutine  ii>  called  by  the 
main  program,  RVR...,  and  is  supplied  with  two 
tracks  as  arguments.  It  does  all  computation 
for  the  spectral  analysis  of  the  filtered  track 
data  and  all  associated  output  (covariances,  raw 
and  smoothed  spectra). 

6.  RVR . . ■  -  The  main  program,  RVR...,  calls  INPUT, 

OUTPUT,  and  FILDEC  to  prepare  the  filtered  track 
data.  It  then  reads  SELECTION  cards  and  calls 
COQUAD  to  perform  the  required  spectral  analysis. 
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B.  Internal  Data  Organization  ' 

The  following  la  a  list  of  the  various  blocks  of 
labeled  COMMON  and  their  associated  variables.  The  table 
which  follows  Indicates  (via  an  X)  which  blocks  are  avail¬ 
able  to  each  routine. 


Block  Name 
DLABEL 
NDMT 

D(20)  - 

RAW 

SMOOTH  - 
NSEL 


Number  of  filtering  coefficients. 

Filtering  coefficients. 

Switch  (true  or  false)  to  control 
printing  of  raw  data. 

Switch  (true  or  false)  to  control 
printing  of  smoothed  data. 

Number  of  SELECTION  cards  to  be 
read. 


DELTAX 

DELTA  -  Distance  between  track  points. 

LAOS  -  Number  of  lags  specified  on  CONTROL 
card . 


LABEL 

NAME(5)  -  COLS.  7-36  of  TITLE  card. 


XXX 

LA(501)  -  Actual  filtering  coefficients  com¬ 
puted  by  CORRET  and  used  by  COQUAD. 


LIST 

IC  -  Switch  (0  or  1)  to  control  printing 

of ^covariance  functions. 

IR  -  Same *as  IC  for  raw  spectrum. 

IS  -  Same  as  IC  for  spectral  estimates. 

I TWO  -  Switch  (true  or  false)  set  by  RVR . . . 
to  tell. 
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QSQ  ■* 

CS(1002)  -  Computational  constants  produced 

QA(501)  -  by  CORRET  for  use  In  COQUAD. 


DLABEL 

DELTAX  . 

LABEL 

XXX 

LIST 

QSQ 

RVR.  .. 

X 

X 

INPUT 

X 

X 

X 

X 

OUTER 

CORRET 

X 

X 

X 

X 

X 

PILDEC 

COQUAD 

X 

X 

X 

X 

X 

X 

C.  Detailed  Description  of  Routines 

Subroutine  OUTER 
ARGUMENTS : 


ARRAY 

-  Floating  Pt.  vector. 

N 

-  Number  of  elements  in 

ARRAY. 

XN 

-  Track  Code:  Track  No 

.  1-A, 

Track  No.  2-B,  • • • , 

Track  No.  S-E. 

KK 

-  1  if  raw  track. 

2  If  smoothed  track. 


OPERATION: 

Outputs  ARRAY,  ten  numbers  per  line,  with  a 
double  space  every  ten  lines  and  a  new  page  every 
fifty  lines.  Each  page  is  headed  by  the  infor¬ 
mation  in  Columns  7-37  of  the  TITLE  card  and  the 
appropriate  track' code. 


Subroutine  CORRET 

PARAMETERS:  (Passed  via  labeled  COMMON) 

LAGS  -  Number  of  lags . 

NP  -  Number  of  filtering  coefficients. 


B-ll  -r 


FACTOR (NF)  -  Filtering  coefficients. 


OPERATION: 

Computes  the  following: 
For  1  •  1,2,  «••,  NF 


NF-i+1 

X 

J-l 


FACTORj  X  FACTOR 


For  1  -  1,2,  ••*,  2  X  LAOS 


0S1  -  cos  C<1-1)  »nm) 


N7 

LA,  •  QA,  +  2  X  I  QA. 
A  l  1*2  1 


For  1  -  2,3,  • * • ,  LAGS  ♦  1 


NF 

LA1  *  QA1  +  2  X  I  PAj  X  Modulo 

J  *  2  X  LAGS]  ♦  1 


Subroutine  FILDEC: 

ARGUMENTS: 

A  -  ARRAY  containing  raw  track  data. 

NA  -  Number  of  points  of  raw  track  data. 

NB  -  Set  by  FILDEC,  number  of  points  In 
smoothed  track  data. 
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OPERATION: 

Smoothes  the  track  data  In  A  and  returns 
smoothed  track  data  to  A. . 


PARAMETERS:  (Via  labeled  COMMON) 

P  -  ARRAY  of  filtering  coefficients. 
NF  -  Number  of  filtering  coefficients. 

For  1  «  1,2,  •**,  NA-NF  +  1 
NF 

Ai  *  Pk  X  Ak+i-l 

NB  is  set  to  NA-NF  +  1 


.  Subroutine  INPUT: 

ARGUMENTS: 

A  -  ,  A  five  column  ARRAY  with  1050  rows 

into  which  track  data  will  be  read. 

I  -  INPUT  will  set  this  to  number  of  points 
per  track  in  the  data  read. 

J  -  INPUT  will  set  this  to  the  number  of 
tracks  read. 


Main  Program  RVR. . . 

VARIABLES : 

A  -  ARRAY  of. track  data  passed  from  INPUT 
routine . 

NDATA  -  Number  of  points  per  track. 

J  -  Number  of  tracks. 


OPERATION: 


See  flow  chart 


Subroutine  COQUAD ; 

Thiss wad  originally  a  routine  for  computing  spectra, 
co-spectra  and  other  related  quantities.  It  was 
obtained  from  the  University  of  California,  Berkely, 
California.  Its  Identification  is  02  BC  COQD.  The 
program  was  written  by  Steward  W.  Smith,  California 
Institute  of  Technology  Selsmologlcal  Laboratory,  and 
was  modified  for  BC  Computer  Center  by  Emily  Harris 
in  February  1963. 

A  flow  chart  for  this  program  is  not  available 
and  operating  Instructions  are  contained  largely  In 
the  above.  We  present  those  quantities  calculated  by 
the  program  which  are  used  by  us. 

The  following  two  terms  are  computed  for  p-0 


to  ffl. 

N-p 

1 

N-p 

N-p 

(TERM  l)p 

•  £ 

1*1 

*1  yl+p 

“  SPp 

£ 

1-1 

*i+p 

X 

1-1 

X1 

N-p 

N-p 

N-p 

(TERM  2)p 

-  £ 

1-1 

yl  *l+p 

•  fpp 

£ 

1-1 

*i*p 

X 

1-1 

y i 

The  following  covariance  estimates  are  computed 
for  p-0  to  m. 

(QX)p  ■  X  formula  for  TERM  1  with  x  replacing 

y  so  that  products  Involve  x*s  only. 

(QX)p  ■  X  formula  for  TERM  1  with  y  replacing 

x  so  that  products  Involve  y's  only. 

(QC)p  "  STTTpT  C(TERM  x)p  +  (TEHM  2V 


(<n)p  *  idbpT £(TERM  “p  2>PJ, 


The  raw  spectral  estimates  are  denoted  (LZ)  where 
Z  may  be  X,  Y,  or  C.  (LQ)  is  computed  by  a  separate 
formula.  These  are  computed  for  h  ■  0  to  a. 

CLI)h  -  i  (9Z)0  ♦  |  “i1  £QZ)p  CO.  i|&  ♦  i  (9Z)n(-l)h 

(M,h  •  i  *jl  <«>p  *ln  ^ 


The  raw  spectral  values  are  subjected  to  a  running 
average  smoothing  using  coefficients  .25,  .50,  .25, 

to  give  the  smoothed  spectral  values  denoted  (WX)h, 
(WY)h,  <WC)h,  and  (VQ>h.  ft 


OPfcJliATI  ON: 


-Ji-15- 


-  *-i‘  i  • 


ENTRY 


•  :Tj 


- Yes/  \ 

STOP.  W— <  NONE=2 


NONE  is  init¬ 
ialized  to  1 
at  load  time 


NONHV  Yes 


CONTROL^ 
Card  / 


NQ 

J 

LAGS  - 
NFMT  - 


NDMT  - 
DELTA  - 

IRAW  - 

ISMO  - 

NSEL  - 


Alternate  Input  Tape  No. 

No.  of  Tracks 
No.  of  Lags 

Switch  (0  or  1)  to  cause 
reading  of  New  Track  Data 
FORMAT  Card. 

No.  of  Filtering  Coefficients 
Track  Data  Point  Interval 
Length. 

Switch  (0  or  1)  to  Control 
Print  Out  of  Raw  Data. 

Switch  (0  or  1}  to  Control 
Print  Out  of  Smoothed  Data. 
No.  of  SELECTION  Cards 


'was  It\.  No 
CONTROL,— « 
v  Card  x 


OUTPUT 
Error 
Message  j 


Road  Filt.  j 
Coef.  Fmt ./ 
Card  / 


Road  / 

Filt.  / 

Cocf.  /  D(l),  .  .  .  .  ,  D(WDMT) 


>4 

■  i 


-  i  i  —  1 7  - 


i  k, 

t .  4 

M 


i 


| 

|  ,i 

t, 
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No  \  Output 
Error 


Call 

IT  2 

COQUAD  with 

h 

— 

Tracks  IT1  & 

IT1 

IT2 

. 

Kl- 

1 

Cal; 

COQUAI 

)  with 

'XIth  1 
- -J 

’rack 

ITWO 

ski  =»  j 

Vi 

. TRUE. 

XJ 

Kl 

IONSEL 


Yes 


o 


>AM 


C  TWO  TPACK  POWEP  SPECTPAL  OeNSfTY 

oexecute  ibjob 

•IP JOB  MAP 

•TB2TC  pvn.n  DECK »2UL 1ST* PEP* 00 
LOGICAL  PAW*  SMOOTH*  fTWO 
DIMENSION  A (1090*9) •  AM(9| 

COMMON  /OLABEL/  NQMT*  0(20)*  PAW*  SMOOTH*  PSEL 
1  fTWO 


LfNOOOOO 
LINOOOIO 
•L|NOOO*0 
A.IST/  fC*  fP*  fS*UNOOOSO 
L1N00040 


DATA  <  AMIfl*  t  •  I*.  J  )  /  IMA* 
1ALP  /  4HSSEL*  3HALS*  3 HALO  / 
tOO  CALL  INPUT t  A*  t«  J  ) 

mdata  •  1  . 

00  200  *  •  1*  J 

12<  paw  i  call  Output <  AM««t* 
IP<  NOMT  «LE*  0  »  60  TO  200 

CALL  2TL0ECI  AII«K|*  Y*  MOATA  » 
200  |2I  SMOOTH  *ANO*  NDMT  «6T,  0  t 


1 MB «  |HC*  I HO*  INK  /•  SCI*  ALS*  LfNOOOBO 

LINOOOSO 

LfNOOOTO 

LfNOOOBO 

LINOOOPO 

f*  ANtKf*  |  |  LtNOOtOO 

LtNOOtIO 

LtNOOtSO 

CALL  OUTPUT t  Atl«K)«  NOATA*  AN<*C* «LtNO0t30 


*2  ) 

LINO0140 

LINO0150 

OO  30«  K  •  1  *  NSEL 

LtNOOtSO 

PEAO  <5*1  )  Tt «  T2*  |C*  IP*  IS*  fTf  *  IT2 

LINOOITO 

IP!  TI  ,NE*  SEL  )  60  TO  900 

LINOOISO 

LtNOOlPO 

IP!  T2  *20*  ALS  I  60  TO  310 

LtNOOEOO 

Linooeio 

IP(  T2  ,20*  ALP  )  60  TO  3*0 

LIN00220 

fTWO  ■  *TPU2, 

L1N00E30 

12 <  IT2  ,6T.  0  )  60  TO  309 

LIN00240 

I TWO  •  *2ALS2* 

LIN00290 

tT2  *  IT1 

L1N002S0 

305 

CALL  COOUAOI  NOATA*  AM*  ITU*  AlftlTSl*  AN<  JTI 1 «  AMHT21  » 

LIN002T0 

GO  TO  300 

LINOOEBO 

*  * 

LINOOE90 

310 

fTWO  •  ,2ALS2, 

LINOO3O0 

00  311  21  •  1 «  J 

• 

LINO031O 

31  1 

CALL  COOUAOI  NOATA*  AI1«K11*  AI1«K1)«  AN<2f1 « 

AH(*t  ) 

i 

LIN00320 

60  TO  300 

LIN00330 

LINO0340 

320 

J1  •  J  •  1 

L1NO0390 

I TWO  •  , TOU2  *  - 

LINO03SO 

OO  321  21  •  1*  J1 

LIN003T0 

22  •  21  ♦  1 

L1NOC3BO 

12 <  J  *N2,  J  )  60  TO  32* 

LIN00390 

22  ■  1 

LIN00400 

I TWO  «  .2ALSE* 

LIN00410 

322 

00  321  23  •  22*  J  . 

- 

L1N00420 

321 

CALL  COOUAOI  NOATA*  All «*I1«  A 1 1*231*  ANlKt ) * 

*Mtr3| 

i 

LIN00430 

300 

CONTINUE 

LINOO440 

GO  TO  100 

LINO0450 

LINO04SO 

900 

WOtTf  (6*2)  2 

LIN004T0 

GO  TO  100 

LIN00480 

1 

20P*AT<AA«A3«3T1 *2161 

LIN00490 

2 

20PNATI23M10N  ATTEMPT  TO  PEAO  TNEI4«T9MTM  MEL 

capo* 

NO  *SEL 

CAPO  LIN0050O 
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END 

•isetc  input  o€ctf*ruLtST«pffp«eo 

SUBROUTINE  fNPUTI  A*  |«  J  ) 

REAL  LA 

LOGICAL  PA W«  SMOOTH 

DIMENSION  All  030*3)  *  EMT(  12)  «  BE«T(t2) 

COMMON  /DLABEL/  NO NT*  0(20) *  RAW*  SMOOTH*  NSEL  /DELTA*/  DELTA* 

1  LAOS  A.AAEL/  NAME(S)  /  XXX  /  LA(SOl) 

DATA  ASTPJC.  DOLLAR*  SLASH*. NOW  IMS*  lH/«  1  /« 

I  EMT(l)  /  |3M(E10,0,A1XA1 )  /«  ILfN*  tP|N  /  AMRLINER*  AH4EINIS  / 

1  XU  NOW  .20*  2  )  60  TO  A 

rr<  now  «eo*  o  )  eo  to  a 
fCTL  •  I 

II  REAO  <S*IOO)  ITST*  NO*  J%  LAOS*  NT  NT*  NOW*  DELTA*  |RAW«  |SHO* 

1  NSEL 

IE(  ITST  »E0«  1EIN  )  00  TO  A 
|E(  ITST  *W,  ILIN  )  60  TO  AO 

PAW  a  .TRUE* 

SMOOTH  a  #TRUE* 

IE(  IOAW  ,20.  0  )  RAW  a'  »EALSE* 

I E(  !SNO  ,20,  0  )  SNOOTH  a  ,2AL*2* 

IE(  NONT  ,L2.  O  »  60  TO  2 

READ  (3*101)  OEMT 

READ  ( 3  *  DENT )  (  D(M) *  Mat*  MONT  ) 

CALL  COROET 

2  trt  NEMT  *6E,  1  )  REAO  (9*101)  ENT 

3  NOW  a  O 

A  REAO  (9*1 09 )  EEY 

A)  REAO  (NO* ) 02)  NAME,  KEYT 

|E (  KEY  ,Nf«  XEVT  )  eo  TO  41 
I  *  t 

3  READ  (NO*EMT)  (  A(|*N)«  N  a  |«  J  ),  CODE 

I  a  I  A  1 

!E(  CODE  ,W.  DOLLAR  *ANO«  CODE  *W.  ASTRIC  *ANO.  CODE  *NE*  SLASH 
1 )  60  TO  3 

I  a  T  -  2 

IE  (  CODE  ,C0,  ASTRIC  )  NOW  a  2 
IE (  CODE  *E0»  SLASH  )  NONE  a  1 

WRITE  ( A* 1 03 )  NO,  J,  DELTA*  I*  NSEL,  LAOS*  NONT 

IE(  NOMT  *LE  *  0  )  60  TO  31 

WRITE  ( A  *  1 0A )  (  0(H)*  Mat*  N©MT  ) 

LAOS  1  a  LA6S  At 

WRITE  ( A  *  1 0? )  (  LA (M) *  M  .  2t  LAOS)  ) 

31  IE(  NO  *  J  *  LACS  *  NSEL  *  IEIX(  DELTA  A  *R  )  *fO,  f)  )  WRITE  (A*1 


LTNO0S20 

LIN00330 

L1N00340 

LINO 0390 

LIN003AO 

LfN003TO 

LIN00380 

LINO0390 

LlNOOAOO 

LINOOAIO 

LINO0A20 

LINOOA30 

L1N00A40 

L1N00A3O 

LINOOAAO 

LINO0AA3 

LINO0A7O 

LINOOAT) 

LIN00A72 

LIN00AT3 

L| NOOATS 

LIN00S7R 

LINOOAAO 

L1N00AR0 

LIN00700 

LINOOTIO 

LIN00T20 

LINOOT30 

LINO0T4O 

LIN00T30 

LIN007A0 

LIN00770 

LINOO790 

LINOOTRO 

LlNOOAOO 

LINOOAIO 

LINO0820 

LIN00830 

LIN00840 

LIN00830 

LINOOAAO 

LINO0870 

LINOOAAO 

LtNOOARO 

LINOOROO 

LINOOOJO 

LIN00O20 

LINO0R3O 

LINOORAO 

LIN00O30 

LINOORAO 

LIN00R70 

LINOORAO 
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109) 

pptupn 

•PITS 

STOP 


<6«10«) 


•pits  <««top) 
1  ! 5 WO*  NSPL 
TCTL  ■  0 
80  TO  11 


1CTL*  1TST •  NQ«  J,  LAOS*  NPNT*  NO«T,  OPLTA*  tPAV* 


Linooppo 

LtNOl OOO 
LtNOIOtO 
LIN0102O 
L1N01030 
L1N01040 
LtNOI04| 
L1N01042 
L tNO 10A3 
LINO I  044 
LtNOl 049 
LtNOl 090 

80PNAT( A8* 12* t 1 *13*11 • f 2*P6*0*2t t • t2t 

80PNAT<12A8t  L1N01070 

80PNAT(6X* 948*38X18)  LtNOl 080 

P0PNA?(49Wl INPUT  PAPANPTPPS  POP  TMP  POLLOWfNG  CONPUTATtONS«/t9MOtLLtN01 090 
IPVATtON  0ATA//3X1 2ttL0Gf CAL  TAPP13//3X0MT9ACKS7Xt2//3X8MDPl.TA  P12*8LtN01 t ©0 
2//3X10MN0.  POINTS 19/21 M0NUM9PP  Op  SSPL  CAPOS 14/I9M0NUP9PP  OP  LAGS t LtNOl 110 
34/26MONUNBPP  OP  PtLTPPtNG  COPP.tA)  L1N01120 

POONATt/lPMOPtLTPPING  COPP*  —  3P17.8//1 19X3P17.8/1X) >  LtNOl 1 30 

F0P9AT( //98MOPPPOP  fN  ABOVP  tNP UT»  ATTPNPT  AT  PPOCPSStNQ  PILL  BP  NL1N01140 
1*08.)  -LtNOl 190 

*0PMaTMh1/1h3//99X.10mTHPSP  TNPN//9BX»3mAPP//93X* I4MTMP  ANSVPPS**LtN01 180 
1 *//90X« 20MN0P  ALL  THAT  PPNA 1NS//99X. 1  ONTO  PP  ©ONp//p9X«2MtS//99X*  L1N011T0 
29MT0  PPCALL //93X *  1 4MTMP  0UPSTt0NS.//89X«8MPVP...  )  LtNOl 180 

P0PNAT1 //29X10M  ACT.  PtLTpP  VALUPS//(Pt9.8*3PtT.8/f X)t  LtNOl IPO 

FOPNATJ 72X18)  LtNOl 200 

POPNATf 1 1 «50M*  SL 1NPP • . CAPO  PXPPCTPD  BUT  NOT  POUND*  CAPO  AS  PPAO/lHL tNO 1 201 

LtNOl 202 
LtNOl 2) 0 
LtNOl 390 
LINOIAOO 
LINO 14)0 
LtNOl 420 
L1N01430 
LtNOl 440 
LINO 1490 
LIN0 1480 
LINO 1470 
LtNOl 480 
LINO  1400 
L1N01900 
LtNOl 910 
LtNOl 920 
L  tNO 1 930 
LINO  1940 
HN01990 
LINO 1980 
L1N019T0 
LINO  1980 
LtNOl 9P0 
LtNOl 800 
L1NO1810 
LtNOl 820 


80 


100 

101 

102 

103 


104 

ITS 

108 


107 

108 
109 

1 08XA8*  t2*tl«t3*2tl *88*0*21 1 *12) 

8NO 

ST98TC  COPPFT  0PCK«8ULtST*PPP«00 
SU8POUT1NP  COPPCT 
OPAL  LA 

CONN  ON  /XXX/  LA(9011/DLA8PL/NP«PACT0Pt23)/T)Pt.TAX/NULL*LAeS 
CONNON  /OSO /  CSt 1002) «0A<901 ) 

LAGSt«LAGS*t 
NPl «NP4 1 
DO  91  1-1 «NP 
0AM)«0*0 
N.NPt-t 
NJ.t 

DO  91  J«1«N 

OAM )«0A < t IPPACTOPt J)4PACT0P«NJ) 

91  NJ«N J* 1 

OAtaOAM  1*0*9 
DO  92  1*1 iLAGSt 

*2  LA  M )«0A  1 

ANG*3* 141 942894 /PLOATtL AOS  ) 

LA6TX0«LAeS4LAGS 
DO  93  1 »1 iLAGTwO 
G>8L0ATtt-1 )4ANG 
93  CStn«COS<G) 

DO  94  t»2*NP 
44  LAM)  •  LAM  )*0A<  t  ) 
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wo  99  Ja2,LA<J31 
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APPENDIX  C 


Two  Dimensional  Power  Spectral  Density  Program 


1.  .  Field  Recording  of  Data 

The  field  recording  of  survey  data  should  follow  a 
pattern  that  simplifies  transfer  of  the  data  to  IBM  cards. 
The  following  pattern  is  recommended. 

The  survey  books  should  have, in  addition  to  ruled 
horizontal  lines, six  vertical  columns.  The  left  column 
is  for  the  x-coordinafce  numbers.  The  regaining  columns 
are  for  elevation  data  in  the  order  below. 


00 

H01 

H02 

H03 

-  Hgij 

10 

H11 

H12 

h13 

Hl4 

20 

H21 

H22 

H23 

H24 

etc. 


On  the  next  pages  the  first  five  lines  of  ground  height 
should  be  completed.  The  x-axls  heights  then  may  be 
read  in  the  left  columns  of  the  first  several  pages.' 

The  lines  of  y-coordinates  5  through  9  should  be  entered 
in  the  five  columns  of  the  next  pages,  etc. 

Suppose  the  coordinate  system  on  the. ground  were 
right  hand  one,  i.e.,  standing  at  the  (0,0)  corner  facing 
the  y-axls  side  of  the  square,  the  x-axls  side  is  to  the 
right.  The  program  prints  out  the  raw  data  and  the 
smoothed  data  in  a  right  hand  system  also,  but  the  y-axis 
heights  are  horizontally  listed.  The  identification  of 
the  ground  directions  with  the  print -cut  is  easily  made 
with  this  in  mind.  The  directional  aspects  of  both  the 
lagged  products  and  spectral  estimates  outputs  follow 
this  pattern.  They  are  in  the  same  sense  as  the  ground, 
but  the  positive  y-axis  is  to  the  right,  the  positive 
x-axls  is  down. 


-'T 


■r*m 
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2.  Computer  and  Operating  System 

The  program  uses  the  IBM  7090  or  709^  computer. 

The  standard  IBSYS/IBJOB  (version  12)  operating  system 
will  operate  the  source  deck  if  the  ‘TIMS'  subroutine 
is  modified  or  dummied  (see  flow  chart).  Since  dif¬ 
ferent  systems  use  different  logical  tape  numbers,  the 
standard  input  tape,  5>  and  output  tape,  6,  may  be  incon¬ 
sistent.  (See  flow  chart  for  appropriate  modifications.) 


3.  General  Purpose 

(A)  Processing  of  a  raw  elevation  matrix. 

Given  as  input  a  raw  elevation  matrix  (adjus¬ 
ted  for  instrument  height)  the  program  computes 
the  following  matrices : 

(a)  Smoothed  elevation  matrix 

(b)  Mean  lagged  products 

(c)  Raw  spectrum 

(d)  Smoothed  spectrum 

(e)  Lagged  products  of  smoothing  coefficients 

(f)  Fourier  transforms  of  L.P.S.C. 

(g)  Corrected  spectrum 

■  Any  of  these  items  may  be  written,  optionally, 

on  SYSOU1.  In  addition,  items  (a),  (f),  and  (g) 
above  are  written  on  an  alternate  output  tape. 


(B)  Processing  a  smoothed  elevation  matrix. 

Using  the  alternate  output  tape  produced  by 
part  (A),  as  input*,  a  subsequent  run  may  be 
made  which  will  produce  items  (b) ,  (c),  (d), 
and  (g)  for  a  selected  submatrix  of  the  smoothed 
matrix,  with  (g)  for  this  submatrix  saved  on 
an  alternate  output  tape. 


•It  is  recommended  that  this  facet  of  the  PSD  program  be 
used  only  in  conjunction  with  input  tapes  produced  as 
alternate  output  tapes  in  part  (A). 


4.  Input  Data  Format 


In  explanation  of  card  preparation  for  all  runs,  the 
following  terminology  will  be  used. 

All  data  will  be  positioned  on  a  card  by  specifying 
column  limits.  'COLS  N  -  M'  means  that  the  data  item 
involved  must  be  punched  in  columns  N  through  M,  N  and 
M  included.  : 

Numerical  data  will  be  of  the  type  floating  point 
(denoted  by  [F])  or  integer  (denoted  by  [I]).  Floating 
point  numbers  must  have  the  decimal  point  punched  and 
lying  within  the  column  limits.  Integer  numbers  must 
be  right  adjusted  in  the  allotted  columns,  i.e.,  the 
units  digit  punched  in  column  M;  the  decimal  point  is 
not  punched. 


5.  Input  Specifications 


Due  to  the  potentially  large  amount  of  data  Involved, 
both  in  size  of  input  matrices  and  number  of  input  matrices 
a  flexible  and,  hopefully,  simple  set  of  input  options 
has  been  provided. 


A  raw  elevation  matrix  Is  punched  in  the  following 
manner,  five  numbers  per  card  in  (F)  format.  The  ele¬ 
ments  of  the  first  row  and  the  first  five  columns  are 
punched  in  the  first  card— 


COLS 


1-10 

element 

of 

row 

1, 

col 

1 

11  -  20 

element 

of 

row 

1, 

col 

2 

21  -  30 

element 

of 

row 

1, 

col 

3 

31  -  40 

element 

of 

row 

1, 

col 

4 

41  -  50 

element 

of 

row 

1, 

col 

5 

In  the  same  manner,  the  elements  of  the  first  five  columns 
and  the  second  row  are  punched  in  the  second  card.  This 
process  is  continued  until  the  elements  of  the  first  five 
columns  and  the  last  row  are  punched.  The  second  five 
columns  are  punched  In  exactly  the  same  manner.  Suppose 
the  number  of  columns  is  not  a  multiple  of  five.  In 
this  case,  the  last  few  columns  are  punched  in  exactly 
the  same  manner  as  above;  e.g.  if  the  number  of  columns 
is  23,  the  last  3  columns  are  punched  using  only  columns 
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1  -  10,  11  -  20,  and  21  -  30. 

/ 

Each  Input  matrix  has  associated  with  it  a  matrix 
title  card  which  contains  a  matrix  number  and  identifi¬ 
cation  information.  The  format  for  this  card  is: 

COLS  1-72  any  alphabetic  text 

73  -  80  matrix  number  (I) 

The  alphabetic  text  in  columns  1-72  will  be  printed 
at  the  top  of  each  output  page. 

Each  elevation  matrix  must  be  followed  by  a  matrix 
termination  card.  The  format  for  this  card  is 

COLS  1-71  not  read 

72  -  72  termination  character 

73-80  not  read 

Each  matrix  which  is  to  be  processed  requires  a 
matrix  control  card  which  controls  tape  assignments, 
output  options,  etc.  The  format  for  this  card  is 

COLS  1  -  6  $GRIDP 

7  -  8  -blank  or  0,  all  input  is  on  the 

system  input  tape. 

-n,  matrix  title  card,  elevation 
matrix  and  matrix  termination 
card  are  on  logical  tape  n, 
in  that  order. 

9-10  «ra,  logical  tape  m  is  alternate 
output  tape. 

11  -  12  -0  or  blank,  this  is  a  raw  eleva¬ 

tion  matrix. 

-1,  this  is  a  smoothed  elevation 
matrix. 

13  -  14  -0  or  blank,  continuous  processing 

of  all  input. 

■1,  after  each  matrix  is  processed 
the  program  pauses  after  print' 
ing  on-line  instructions  to 
the  operator  to  allow  continu¬ 
ing  or  terminating  the  run. 
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15  -  16  -0  or  blank,  for  all  matrices  In 

all  runs,  with  the  following 
exception:  Suppose  3  matrices 
have  been  processed  and  their 
output  written  on  an  alternate 
output  tape,-  The  current  run 
Is  to  process  several  more 
matrices  and  It  Is  desirable 
to  write  the  alternate  output 
from  this  run  on  the  same  tape 
reel,  after  the  output  from 
the  first  3*  In  this  situa¬ 
tion,  this  field  on  the  $GRIDP 
card  for  the  first  matrix  only 
Is  punched  01. 

17  -  18  or  blank,  the  elapsed  time  is 

printed  after  each  of  the  items 
(a)  through  (g)  is  computed. 

*0  or  blank,  no  timing  Is  done. 


The  following  fields  control  the  suppression 
of  output  on  SYS0U1  of  the  Indicated  items. 

0  or  blank  suppresses  output  of  that  particular 
item  while  01  causes  the  item  to  be  written. 


COLS  19 
21 
23 
25 
27 
29 

31 

33 

35 

73 


20  elevation  matrix  as  read. 

22  smoothed  elevation  matrix. 

24  mean  lagged  products. 

26  raw  spectrum. 

28  smoothed  spectrum. 

30  lagged  products  of  smoothing 

coefficients . 

32  Fourier  transforms  of  L.P.S.C. 

34  corrected  spectrum. 

72  not  read 

80  matrix  number  (I).  The  program 

will  search  the  tape  which  is 
specified  as  being  the  matrix 
tape  (if  cols  7-8  are  0  or 
blank,  SYSIN1,  if  cols  7-8-n, 
logical  tape  n)  for  a  matrix 
title  card  with  this  number 
in  cols  73-80. 


re** 
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Data  Interval  Card 

COLS  1-10  column  Interval  (F) 

11  -  20  row  Interval  (F) 


Matrix  Size  Card 

COLS  1  -  10  number  of  columns  in  the  eleva¬ 
tion  matrix.  (I) 

11-20  number  of  rows  in  the  elevation 
matrix.  (I) 

21  -  30  value  of  p.  i.e.  (the  number 

of  columns  in  the  smoothing  co¬ 
efficient  matrix  -  l)/2.  (I) 

31  -  40  value  of  o.  i.e.  (the  number 
of  rows  in  the  smoothing  co¬ 
efficient  matrix  -  l)/2.  (I) 

41  -  50  value  of  MX.  (I) 

51  -  60  value  of  MY.  (I) 

The  corrected  spectrum  will  be 
(2  •  MX  +  1)  by  (2  •  MY  +  1) 


The  Q-  matrix  or  spectral  smoothing  matrix  (3  by  3) 
la  presented  on  punched  cards  in  the  following  format 

COLS  1-10  element  in  row  1,  column  1  (F) 

11  -  20  element  in  row  1,  column  2 


61  -  70  element  in  row  3*  column  1 
second  card 

COLS  1  -  10  element  in  row  3,  column  2 

11  -  20  element  in  row  3,  column  3 


The  B-matrix  or  matrix  of  smoothing  coefficients 
has  2p  +  1  columns  and  2a  +  1  rows.  As  it  is  symmet¬ 
rical  about  the  origin,  only  the  upper  half  of  the 
plane  >  0  is  specified  in  (2o  +  1)  *  (p  +  1)  numbers. 
That  is,  columns  p,  p  +  1,  .  .  .  2p  +  1. 


Punching  by  columns,  7,  10  column  (P)  fields  per 
card,  column  p  is  punched,  then  column  p  ♦  1,  etc. 


Finish  Card 

The  format  for  this  card  is: 
COLS  1-6  $FINIS 


6.  Peck  Set-up  -  Single  Elevation  Matrix 

With  the  above  mentioned  cards,  a  complete  deck, 
set  up  to  process  one  raw  elevation  matrix  consists  of 
the  following  cards  in  the  Indicated  order.  This  Ulus 
trates  deck  set-up  If  no  alternate  input  tape  is  speci¬ 
fied. 

$JOS 

Installation  identification  card. 

$ EXECUTE  IBJOB 

$IBJ0B 

BINARY  OR  SOURCE  DECK  FOR  PSGRID  PROGRAM 


$DATA 

$GRIDP 

Matrix  title  card  * 

Data  interval  card  f 

Matrix  size  card  # 

G-raatrix  # 

B-raatrlx  § 

Elevation  matrix  » 


Matrix  termination  card 
$FINIS 


If  an  alternate  Input  tape  is  specified,  the 
*'ed  items  do  not  appear  on  the  system  input  tape 
but  are  placed  on  the  alternate  unit  in  the  follow¬ 
ing  order. 
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Matrix  title  card 
Elevation  matrix 
Matrix  termination  card 


7«  Deck  Set-up  -  More  Than  One  Elevation  Matrix 

To  process  more  than  one  elevation  matrix  per  run, 
the  termination  characters  on  the  matrix  termination  cards 
control  the  Input  for  the  next  matrix  to  be  processed. 

•  In  column  72  -  This  matrix  Is  the  last  matrix  to 

be  processed  in  this  run.  Execu¬ 
tion  will  be  terminated  after 
this  matrix  is  processed. 

/  in  column  72  -  In  processing  the  next  matrix, 

fed  items  are  to  be  read  In  again. 

$  In  column  72  -  In  processing  the  next  matrix, 

I'ed  items  will  be  the  same  as 
for  this  matrix,  therefore  these 
items  are  not  included  after  the 
IQRIDP  card  fcr  the  following 
matrix . 


The  last  card  of  every  data  deck  must  be  a  $FINIS 
card.  There  should  be  only  one  $FINIS  card  for  each  run. 


Output  Specifications 


Alternate 


The  information  written  on  the  alternate  output  tape 
during  t*  processing  of  a  raw  elevation  matrix  is  organ¬ 
ized  In  '  .j  following  manner: 
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Matrix  Title  Card  (as  read) 

Data  Interval  Card  (as  read) 

Matrix  Size  Card  (as  read) 

except 

COLS  1-10  number  of 
columns  In 
smoothed 
matrix 

11  -  20  number  of 
rows  In 
smoothed 
matrix 


G-matrix 

Smoothed  Matrix* 

Pourier  Transforms  of  L.P.S.C.** 

Matrix  Termination  Card  -  Termination  Character/ 

Matrix  Title  Card 

COLS  1-72  (as  read) 

73  -  80  (as  read)  +5  , 

Corrected  Spectrum  Dimensions 

COLS  1-10  number  of  columns  In  correc¬ 
ted  spectrum 

11  -  20  number  of  rows  in  corrected 
spectrum 

Corrected  Spectrum** 

Matrix  Termination  Card  -  Termination  Charac¬ 
ter  $  except  for  last 
matrix  in  run ,  in 
which  case  termina¬ 
tion  character  is  * 
and  COLS  73-80  con¬ 
tain  99999999 


Due  to  the  problems  created  by  the  necessity  of 
finite  representation  of  numbers,  the  three  *'ed  items 
are  written  on  the  alternate  output  tape  in  a  manner 
which  will  preserve  the  entire  machine  representation 
of  their  values.  The  format  for  all  of  these  numbers 
if  +  0 . zxxxxxxx£+yy .  This  notation  represents  the  num¬ 
ber"” 
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♦  0 


.txxxxxxx  X  10  +  yy 


Numbers  written  in  this  notation  always  require  15 
columns  and  z  is  always  non-zero,  unless  the  number 
Itself  is  zero. 

*  The  smoothed  matrix  is  written  in  a  manner  analogous 
to  that  of  the  raw  matrix,  l.e.,  4  columns  at  a  time, 
the  values  occupying  COLS  1-15,  16  -;30,  31  -  45  and 
46  —  60. 

**  The  Pourier  transforms  of  L.P.S.C.  and  the  corrected 
spectrum  are  written  in  a  slightly  different  manner, 
l.e.,  placing  4  numbers  per  card  image  the  entire  first 
row  is  written.  Starting  with  the  next  field,  the  sec¬ 
ond  row  is  written.  This  process  continues  until  all 
rows  have  been  written. 


9.  Deck  Set-up  -  Processing  a  Submatrix  of  a  Smoothed 
Matrix 


Using  the  alternate  output  tape,  described  In  8., 
as  the  alternate  input  tape,  any  submatrix  of  a  smoothed 
matrix  may  be  processed  as  indicated  in  3.(B).  Deck 
set-up  is  as  indicated  in  6.  and  7.  with  the  following 
exceptions: 

The  matrix  control  card  (I0RIDP  card)  must  have 
'01'  punched  in  columns  11-12. 


All  fed  and  *'ed  items  In  6.  are  removed  and 
replaced  by  the  following  submatrix  card. 


COLS  1-10 
11  -  20 
21  -  30 
31  -  40 


lower  column  limit* 
upper  column  limit 
lower  row  limit 
upper  row  limit 


•Both  rows  and  columns  are  numbered  starting  at  1,  even 
though  the  print-out  in  the  raw  matrix  run  begins  num¬ 
bering  at  0. 
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SUPPLEMENT  A 


Two- Dimensional  Power  Spectral  Density  Program  Notes 


The  PSD  program  Is  organized  as  a  main  program  and 
four  subroutines.  Labeled  common  Is  used  where  appro¬ 
priate  and  Identical  variable  names  are  used  when  Items 
are  transmitted  via  labeled  common.  See  Supplement  B 
for  a  list  of  the  Important  variable  names  and  their 
usage. 

The  following  Is  a  brief  description  of  each  of  the 
decks  which  comprise  the  PSD  program.  Plow  charts  are 
attached  for  the  main  program,  deck  name  IN,  and  the 
major  computational  subroutine,  deck  name  FSGRID.  It 
is  assumed  that  the  source  listings  of  the  remaining 
decks  are  self-explanatory. 

'CLOCK'  A  machine  language  program  which  places 

the  integer  number  representing  the  remain¬ 
ing  execution  time  (in  .6  seconds)  in  the 
AC  each  time  it  is  called. 

'TIME*  A  FORTRAN  subroutine  which  prints  the  total 
elapsed  time  (time  since  last  call  with 
zero  argument)  and  the  time  elapsed  since 
the  last  call,  in  minutes  and  seconds. 

NOTE:  Both  of  these  decks  may  be  removed 
and  a  dummy  'TIME*  deck  substituted  if 
the  user  so  desires. 

'WRITER'  A  FORTRAN  subroutine  which  writes  the  argu¬ 
ment  matrix  on  the  output  tape,  6,  5  columns 
at  a  time,  indexing  the  rows  and  columns 
as  specified  by  argument  values. 

'PSGRID*  A  FORTRAN  subroutine  which  computes  and 

calls  'WRITER'  to  write  items  (b)  through 
(g) ,  given  the  smooth  matrix. 

'IN'  The  FORTRAN  main  program  which  does  all 

Input  for  both  raw  and  smoothed  matrix  runs, 

'  computes  the  smoothed  matrix  if  necessary 
and  calls  'PSGRID'  for  the  remaining  com¬ 
putations  . 


rr'“-»-  — "^W"-  #"•**■»•» > '%*’’ ir.-n nr  -- 
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Included  with  the  source  deck  of  the  PSD  program  are 
two  file  specification  decks.  These  assume  that  FORTRAN 
logical  tape  7  will  be  designated  as  the  alternate  input 
tape  and  FORTRAN  logical  tape  8  will  be  the  alternate 
output  tape.  Both  tapes  are  made  up  of  80  character 
logical  and  physical  records,  i.e.  card  images.  This 
structure  must  be  used  to  be  consistent  with  READ  and 
WRITE  statements  throughout  the  PSD  program. 


Single  Location  Variables 
NX  *  Nx  +  1,  the  width  of  the  H  matrix. 

NY  ■  Ny  +1,  the  height  of  the  H,  matrix. 

NRHO  -  p,  the  horizontal  limits  of  the  B  matrix 
NSIG  *  o,  the  vertical  limits  of  the  B  matrix. 
N3R  ■  2p  +  1,  the  width  of  the  B  matrix. 

NBS  »  2o  +  1,  the  height  of  the  B  matrix. 

NHR  -  p  +  1 
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NHS  "o+l 

NSHX  ■  N  +  1.  2p  *  n  +  1,  .  the  width  of  the  h  matrix. 

A  •  A 

NSHY  •  N  +  1  -  2o  »  n  +  1,  the  height  of  the  h  matrix. 

*  I  * 

MX  ■  mv,  the  horizontal  limits  of  the  <j  and  f  matrices. 

'  I 

MY  ■  m  ,  the  vertical  limits  of  the  o  and  f  matrices. 

7  i 

NGX  «  2m  +  lt  the  width  of  the  matrix. 

*  i 

NFX  -  NGX 
FX  *  NGX 

NFY  ■  2my  +  1,  the  height  of  the  f  matrix. 

FY  »  NFY 

NGHX  -  mx  +  1  | 

NGHY  -  my  +  1 
NFBX  -  2mx 

NFBY  -  2m  ! 

y  i 

N2R  »  2  j 

| 

NWX  ■  4  +  1,  the  width  of  the  W  matrix. 

DX  *  x,  the  horizontal  data  interval. 

DY  *  y,  the  vertical  data  interval. 

INTAPE  »  the  alternate  input  tape  numbei . 

OUTAPE  ■  the  alternate  output  tape  number. 

FLAG  »  0  for  raw  run,  /  0  for  smoothed  run. 

IBUG  *  switch  for  on-line  messages. 

IBTAP  ■  switch  for  previously-processed  alternate  out¬ 
put  tape. 

I 

I 

j 
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ITIME  *  switch  for  timing  print-outs. 

IPT(l)  through  IPT(8)  ■  printing  switches  for  raw  ele¬ 
vation  matrix  and  items  (a) 
through  (g)  [listed  under  3*  (A) 
respectively . 


rab 


(n  -  a  )  (n  -  b  ) 
x  y 


n  n 
x  y 

2  2  h,  .h,  .  , 

k=0  j=0  k+a»J+b 


where  n  =  N  -  2p+l 
x  x  r 


n  =  N  -  2ct  +1 

y  y 


and 


-m  <  a  <  m 
x  —  —  x 


-m  <  b  <  m 

y  -  -  y 


m  m 
x  y 


a& 


AxAy  2  2 

a  =  n  b  =  m 


x 
aa 


72b  COs[2,(?HrTr  + 


2  m  TT}  1 

y 


and  -m  <  a<  m 
x  -  -  x 


-m  <j3  <  m 

y  -  -  y 


f 


1 

=  2 


1 


nv  2  =  -  1  j-f-1  J2j  fn+2. 1/H 


and  if  rj  +  1  or  v  +j  occurs  out 
side  the  border,  take  the 
reflected  index  inside. 


C  APFA  SPfCTPAL  DENSITY  PP06PA* 

*e*ffcur»  ipjcp 

*TPJOP  W.*P 

S IflFTC  PSGPTO  0e,Ctr«00»FULTST»PeF 
SU8P0UTTNF  PSGPTO 

C  MOOfFteO  SPffCTPAL  OFNStTY  FPOG  STAPTTO  Irt  4F  T  1*64 


INTCGFP  OUT  APE  «  Ft- AG 


COMW*:  /AOOAYS/  SMft05«10O>* 

1  /St?€S/  *X«  "*•  M*W* 

2  /CONTBL/  T*T<ai«  FLAG* 

3  •  /LAP/  LASFL<!2»«  N|0 


PHt<30«231»  P(23«2P»«  G(3«31 
NBS*  NSM*«  NSMV*  *FX«  n*»  DY«  NOHO 
1 T  f FF *  OUTAPF 


OIMFNStON  GA(S0«23) «  F(30«50)«  FtJ<SO»SO»*  t#<30«23),  FSTS0.50). 
t  hFA0<6«6> 


DATA 

1 

.2 

3 

A 

S 


H€A0<1 ill 
H3A0M *2) 
HFAOIl *3) 
HFAOtl «4) 
HFAO<1 *5) 
HfAOI1.6) 


/  34H«FAN  LAGGFP  PWOfHJCTS*  (COVAOl 
/  34MOA*  SPFCTPUM.  <F1 
/  34HSPOOTHFO  SPFCTOUP*  IF-PAOT 
/  34MLAGGFO  PPOOUCTS  OF  SMOOTHING  CO€F • 
/  34HF0UPTFP  TPANSFOPPS  OF  L«P«S*C« 

✓  34MCOPPFCTEO  SPFCTPUP*  fF-STAP) 


/♦ 

/« 

/« 

/♦ 

/« 

/ 


EOUtVALFNCF  (FS.FP1«  <GA»W1 

Cl 

c  PFAO  IN  INPUT  DATA 

PI  2  *  2«0*3»t41Fo 

c  COMPUTF  ThF  mean  LAGGFO  PPOOUCTS 

c  A  LESS  than  7FPO 

212  NGHX.MX+1 
NGMY«PYa1 
NGX*NOHX*PX 
D0127I8*t  *PX 
KSO.MX+l-ta 
X67-NSHX-F6P 
OOl 27J8* 1 .NGHY 
K70*J8-1 
K68*NSHY-X70 


(COVAPTFNGFSl 


C13A 

SUH*0«0 
0012816*1 *K67 

|7aKftO+l«. 

O0l28J7*l «X6P 
J8*K704J7  • 

TFMP«SMI I7,J7)*GH( 16* J61 
128  SUN-SUNtTFMP 

0*K67*K6P 

127  GAII8*J8) * SUM/O 

C14  A  FOUAL  TO  AND  GOFATFP  THAN  7800 
213  001T7I8.NGHX.NGX 
X6P.I8-PX-1 
K7t«NSMY-<6« 


GPOOOOOO 

goooooio 
GP000020 
6P000030  ; 

60000040  f 
G0000030  i- 
60000060 
00000070 
GPOOOO80 
60000000  ^ 
GOOOOIOO  ■; 
GP0001 10  i 
GP000120 
60000130 
GR000140 
GOOOO ISO 
GOOOO 1 flO 
GP000170 
GOOOO ISO 
GPOOOIPO 
GPOO0200  ? 
GPOO0210  i 
G0000220 
G0000230 
GOnO0240 
G0000230 
GPOO0260 
GPP00270 
GP0O028O  . 
GP0002O0 
GPD0O3OO 
GOOOO 31 0 
60000320 
60000330 
GPOO074O 
60000330 
GOOOO 360 
G0000370 
GOO 00380 
G000070F 
60000400 
GOOOOA 1 O  | 

Gonno*?o  5 

GPOO0430  < 
60000440  , 
60000430  | 

. GO0O048O  ’ 
60000470  j 

G000048A  ! 

6Ppo04QO  ; 

ooonosno 
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00177 J8«l  »NGMY 

W70«J8-1 

K72*NSHV-K70 

CI4A 

SUB«0«0 
DC17B17-1 «K71 
I6-I74W60. 

00178J7«1 «K72 
J6«J4*4*70 

TpBP«SH<  17.  J7)*SH(  16*  J6> 

ITS  SUB«SU*4TFMP 
0»K71  *K7? 

177  GAtlB.jai-SUB/O 
CIS 

IF<  IPT(3)  *NF*  o  ) . . 

1  CALL  WPITFP {  fiA*  SO*  O*  BV«  -BX*  BX*  MFAPH*!)  ) 

I F<  IT  1 PF  «NF«  0  1  CALL  TtAfftl 

C16 

C  PAW  SPECTRAL  CALCULATIONS 

215  NFY«2*PY41 
Nry*2*«y+i 
fy«nfY 
fx.nfx 

lppx  .  px  ♦  l 

Lppv  ■  BY  ♦  1 
D0131I10«1«nfx 
alpha  ■110-  LPPy 

ALPHA  ■  ALPHA  /  FX 

00131 J10«l «nfy 
BFTA  ■  J!0  -  LPPY 
BETA  ■  BFTA  /  FY 
SlP»0*0 

C17 

0013619*1 «N6X 
A  ■  IP  -  LPPX 
A  ■  A  •  ALPHA 
ooi32J9»2«nghy 
PF«  J9- 1 

A NO 1  •  PI2  »  I  A  ♦  »  BFTA  | 

C17A 

TpPP«  6A< T9« J91*C0S<AN61 > 

132  SUP-SUP+TSPP*?.0 
AN62  ■  PI 2  *  A 

136  SUN«SUH46Af  T9*n*C0Sf  AN62) 
t  31  F<I 10* J10)»DX*OY*SUP 
Cl  8 

IF  1  fPT(A)  « NF  «  0  >  ...  - 

1  CALL  WPITFRI  F»  BO «  -«Y«  BY ♦  -By,  BX*  HFA0ft«2>  ) 
IF(  ITIBF  »NF,  o  )  CALL  TtBFft) 

CI9 

C  SBOOTHFO  spfctpup 

C  INSIOF  PATPIX 

217  NFHX»NFX-J 
NFBY»NFY-1 


gpooosio 

6PD00S20 
0*000530 
6P000540 
GPOOOSSO 
GRD00*6o 
6*000*70 
6PD00580 
GR000890 
GRD00600 
GR00061O 
•  6POO062O 
6R000630 
6PD00640 
6POO06B0 
GR000660 
GR000670 
6RPOO6B0 
GRD00690 
GP000700 
GP00071 0 
GP000720 
GPD00730 
6R000740 
GPOOOTBO 
G0000760 
GPOO0770 
GR000780 
GR000790 
GROOOBOO 
GPOOOB10 
GR000B20 
GRD00830 
6ROO084O 
GROOOSBO 
GP000860 
6P000870 
6PD00880 
GP00089D 
GR000900 
G»O0091 0 
6*000930 
6PDO0930 
GP000940 
6POO09BO 
GR000960 
CR000970 
6P0009B0 
GR000990 
GRPOl 000 
GPPOIOIO 
GRPOl 020 
GRPOl 030 
GR001040 
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D0t40113*2*NP8X 

00140J13*f*N88Y 

SUWaO.O 

00141111*1*3 

tt2*113*m-2 

00141 J11«1«3 

J12*J!3* J11-2 

cao 

TnPt«(ttiuiittr(fttiJUt 

141  SU«*SU**T88P 
140  2Bt1lS*J131*9UM 

C21 

218  ooi43ii3«**Nrex 
STJW«0»0 
SUWAaO «0 
001431 1 4a t  *3 
M8*113*t14-2 

TE*P*G1  !14*n*2<ltS.21-*G<  t 14*21*41118*1 1441114* 3>*4l  1 18*21 
TEM84*61t14« 1t*41tlS*N48Y14«1t1**21*41t1S«M4Y>4d1tt4*31 
X*rtt 15«NESY> 

suwasuiit4Tnap 

143  SU«A*SUWA*TE88A 
20 1113*1)  *SU8 
1*2  EBtl13*MEY>*9U8A 

C22 

214  00144J1 3*2«NEBY 
SUM*0*6 
3U«A*0*0 
D01*8Jt 441 *3 
J!S*Jl**J13-2 

TEWPaGI 1 «Jl*l*E12.JlS1*<*t2« J1«|*E11 , jt8»*4tS*Jt4|«4t2« Jt8> 
TEMPAaGIl  «  Jl41*T(Nr«X«Jt 914412*^1 41  *4 1Mrx«J13><4413«J141 
X*E(*EBX«Jt8) 

SU**SU84>TEW» 

143  SWA*SUWA4TE8PA 
E811 «J13)*8U8 

144  EB1NEX* JtSlaSUMA 
C23 

220  K27»2 
K28*1 
K24*2 
K30*1 

001 *6C 1 4a 1 *4 

SUW*0*0 

X26**24 

001471.20*1  *3  ...  ..  -  - 

TEMP*0(L20*1  >*E(K26*K27»*«a.f0«2l*41K26*Kf8»4«1t.20*3)*EtK26«K27> 
C24  .... 

IE (L 20-2 1 1 48 » 1 44* 1 47 
14®  K2«*K30 
GOTO 147 
1*9  X2*«X29. 

GOTO 147 

1*7  SU4aSU*4TE*P 


60001080 
60001060 
6000 1070 
60001080 
60001040 
60001100 
60001110 
60001120 
60001130 
60001140 
60001 190 
60001160 
60001170 
60001180 
60001140 
60001206 
60001216 
60001220 
60001230 
60001240 
60001286 
6000 1260 
60001270 
60061280 
60001240 
G0001306 
6000 1 31 0 
6006 1320 
60001330 
600013*0 
68001380 
68001360 
68001370 
G0001380 
GOOOI 340 
68001400 
G0001410 
68001420 
68001430 
68001440 
68001*80 
68001*60 
G8001470 
68001*80 
68001*40 
68001800 
68001810 
68001820 
68001830 
64001840 
68001880 
68001860 
68001870 
68001880 


••• 


C  -25 


GOTO <190* 191 .19*. 194) *L19 

190  fH!<J  «1)«SU* 

»cs9«Nrex 
K30*HTX 
GOTO 1 40 

191  PBJNPX.I )«SU* 

K29«2 

*30*1 

k*7«npby 

k«»mfy 

GOTO 140 

e*9 

19*  roii»Mry>«sm* 

«?***•  NTBX 
K30«NPX 
GOTO 140 

194  28<NFX«NrY)«SUM 
140  COWTtMJE 

err  .  - 

tr<  ipt(9»  .Mr*  oi  _ '  '  '  '  ~ 

1  CALL.  W9TTEP<  FB.  SO.  >*«Y«  *7.  -***•  **•  HEAOM.3) 
TF(  ITIPE  .HE.  0  )■  CALL  TIPEM) 
ir  <FLA#t  222.  m*  22« 

£2g 

c  COMPOTE  CORRECTION  COEFFICIENTS 

C  U  LCSS  THAN  ZERO 

JM  W29»MPM0-*N9H0 
001 STt 1 4*1 »N29 
K38«N2R4l -1 1 8 
K40.NBR-K38 
D01S7J18«1 «NBS 
r30*J18-l 
XAI •NBS-K3P 

CEP 

SWM*0.0 

0019§tt7»1 «X40 
1 1«*1 174*39 
00188J19-1 »**i 
J17«J194*39 

TEPP«S( tt7*J17)*9<T19»Jt9) 
tse  sum«su*ttemp 

197  »U1S« J181-SU* 

C30  _  .  _ 

c  u  equal  to  op  opeatep  than  zepo 

NWY-NCP-.M2P  . 
tz 3  DOISPT 18*N9P»N9Y 

K42* t 1 8-NBP 

K40»NBP-**2 

5  001S9J18-1 .MSS 

K3P»Jie-i 
K41»NBS-*39 


SUM«0.0 


69001990 

69001000 

G9OO101O 

GPOO 1020 

GPOO 103© 

SPOOl 840 

0P001890 

GPO01880 

GP001870 

GPOO180O 

69001090 

OPOOITOO 

GPO01710 

6900172© 

GP001730 

0P001740 

6P001790 

0PO01780 

6PD0I770 

6PO01780 

GPOOITPO 

GP001800 

6P001810 

69001820 

6POO 1 830 

GPO01840 

GPOO1890 

GPOO 1 800 

GPO01870 

GPOO  1 880 

69001890 

GPOO 1 POO 

OPOOIPIO 

GPOO IP* 0 

GPOO 1930 

GPOO 1940 

GPOO 1990 

GPOO 1980 

GPOO 1970 

GPOO 1990 

GPOO 1990 

gpooeooo 

69002010 
GP002020 
69002030 
69002040 
GP002090 
6900*080 
GRD02070 
GP002080 
69002090 
GP002100 
GPO021  tO 
69002120 


C-28 


00160t1?«t«K40 
llOaf 174*4* 

00t60J17«l «*41 
J17aJ174*30 

Tp8Pa8tt*7*Ji7)*8ms*Jt7» 

140  SUMa9U84T*MP 
187  W<118*Jt8)aSUM 

51 

tr(  107161  tNP«  0  1- 
I  CAU.  W*1TEB(  W*  SO*  0*  M8S-1* 

twt  iTise  *ns*  o  »  call  Tfxrm 

52 

LOW*  •  ox  ♦  t 
its  00168120a! «NTX 
S  a  120  -  LP8X 

s  a  s  /  rx- 

OOieSJSOaf  *NGMV 
TaJSO-1 
T  a  T  /  FV 
9UWa0«0 

00107  121a) «NXX 
Ua  121 -MSP 
U  a  U  *  S 
00166- J21a2* MBS 
V.J21-1 


A HOI  a  012-4  (  V  -  V  •  T  « 

TE**»aWt  **COS  tANOI) 

100  SU8aSU84TS«P*2*0 
AMO?  a  Ptt  *  U 

107  SUMaSUt*4WU*1*n*eOS  <4*821 

tr<ASS<SUM)*LT«  0*00001)  SUMaSTQN(e*OOOOOt*SU«) 
10S  PMt  M  20  «  J20  1  aStJM 


HTAOll *«) 


XPtTS<0UTAP2*110t  C  (  PMt  f  1  ♦  J)  *  I  •  t*  "0*  )•  j  •  1*  *6W  1 
•PITS  lOtfTAPS *  1 03) 

220  T Pi  1PT<7)  *M2*  0  1 

1  CALL  WPfTEOf  PMt*  SO*  0*  MOMY-I*  «•*»*«  ) 

f*<  ITIMS  «MS*  O  )  CALL  TfWSMl 

C34A 

C38 

227  00167t22al*M*X 
0010PJ22al  «Miry 
123»T22 

J23aM0MYal-J22 
ir(J22*0T*MY)  J23aJ**-**Y 
100  PS< 122* J22)aFB( 122* J22) /PM» 1123* J*3| 

KSOaMOMXal 

KSlaMOMV*>l 

r*lNOM*«MOMr)atrS»MX«KSl  >4PS<H0MX**S1  »4««90.W1  )4PSMCS0«M0MY))/ 
*4*0 


1*  MOW 


WOITC  «OUTAP€»101) 
woire  <outap€*iio) 


LAB2L*  MfO*  M2*  07 
<  1  FJUtJl,  1  a  1* 


00002*130 
000021-40 
0OO02[tS0 
000021 00 
00002*1 70 
00002180 
00002100 
00002200 
00002210 
00002220 
00002230 
00002240 
GOD022S0 
00002280 
00002270 
00002280 
G0002200 
00002300 
G0002}310 
G0002*320 
00002*330 
80002*340 
00002)380 
00002*380 
00002*370 
00002380 
00002(300 
00002*400 
0000241 0 
00002*420 
00002)430 
00002)440 
00002*480 
GOO0240O 
00002470 
GOO02480 
00002440 
G0002800 
00002810 
00002820 
G0002830 
00002840 
00002880 
GPO02860 
G0002870 
00002880 
G0002840 
00002600 
00002810 
00002620 
00002630 
00002640 
00002680 
G0002660 


1E(  1PT<8)  |W|  O  I 

I  call  wwjTmt  rs«  so«  -*r*  by*  -*rx*  bx«  heao<  i  *«>  t 
tr<  trim  *ne*  o  j  call  tihemi 
return 

101  EOBHATT 12AS* TS/2T10««0XA1) 

ICt  B0RRAT<7B10.3«10X) 

103  EOBHATT 71X1 M/RXA 11 

no  eobhat|4eis.8«2oxi 
eno 

•Torre  jn  deck«oo*xulist»ree 

INTEGER  OUT  APB *  FLAG 

COHHOM  /ARRAYS/  SMT105*100)*  PHtT50«28|*  8!t5*25)«  6»3*3> 

1  /SIZES/  NX*  BY*  NOB*  NOS*  NfWX*  NS«¥*  HEX*  0*«  OV* 

2  /CONTRL/  1 PT <81*  FLAG*  ITfRE*  OUTAPE 

3  /LAO/  LABEL! 121 •  NTO 

OfHENSTON  M T 109* 1 001  *  HEA0(6«31*  ELAOELC 121 
EOUTVALBNCE  TM»SM)«  TTO.NtOI 

DATA  1CTL#  AST «  SLASH*  OOL.  TS*T«  JEC«  tf»lN  /  GMROBTOP*  1 

1  1M«*  1*  |«  0HAP1N1S  /« 

2  MEAD! 1*1)  /  32HTHE  SMOOTHING  COEEETC TENTS  /« 

3  HEA0T1«21  /  32MTHE  ELEVATION  HATBTX  /• 

4  HEAD! 1*3)  /  32HTHE  SMOOTHED  MATRIX  / 


CALL  Tt«E!01 

lOOO  READ  (5*100)  t,  JNTAPE*  OUTAPE*  FLAG*  T8UG*  TOTAP,  IT1ME* 
1P(  1  ,E0«  IPtN  )  GO  TO  500 
TE <  t  «ME.  TCTL  )  GO  TO  200 

IE 1  TOTAP  «E0*  0  «OR»  OUTAPE  *E0,  «  >  GO  TO  1002 
tool  RE AO  (OUTAPE* 10«)  ELAOEL*  COoE 

TE!  COOE  *NE.  AST  )  GO  TO  1001 
BACKSPACE  OUT  A PE 
WRITE  TOUT APE • 306 ) 


1002  JEC  •  1 

TE<  I BUG  *£9«  0  1  GO  TO  1 
PRINT  100 
PAUSE 

CALL  SSWTCHT  «*  n 

TE(  I  ,E0.  1  1  60  TO  501 

1  ,/TAPE  •  5 

TE(  INTAPE  *6T.  0  )  JT APE  ■  1NTAPE 

2  BEAD  ( JTAPE* 101)  LABEL*  TOT 


GROO2670 
GR002000 
GR0020R0 
GRO 02700 
GR00271 0 
6R002720 
GR002730 
GPO 02740 
GRO027B0 
6R002700 
GR002770 
GRO027R0 
GRD02000 
6BD02810 
HRhO  6P002020 
GR002830 
6R002840 
GBO028B0 
6R0028G0 
GRD02870 . 
GRO02880 
-GBD028B0 
H»*  tM/*GR002900 
6R002R10 
GRO02R20 
GRD02B30 
GR002940 
6O002R50 
6B002R60 
5RO02R70 
6RD02R80 
TPT*  TDGR0029R0 
G»003000 
GRO03010 
GROO302O 
GRO03030 
GRO03040 
GR003050 
GR003060 
GP003070 
GP003080 
GR0030R0 
GR003100 
GB0031 10 
GR003120 
60003130 
GR003140 
0RO03150 
6R003180 
GR003170 
GRO03180 
6RO031R0 

60063200 

GR003210 


C-28 


TEC  TOT  *NE*  fO  *  60  TO  S 

(.TAPE  •  9- 

JE<  Ft  A  6  *NE.  0  I.  LTANE  ■  TNTAN1T 

ie(  tsvt  .n«  o  »  so  to  *i 

XEAO  (LT AXE* 108)  ox*  or 

XEAO  (LTAXE.T03)  NX*  NY.NXwO*  MSt«*  NX*  NY 

XEAO  (lTAK*|Ot)  C  <  «(t*JI*  I  •  I*  3  >«  J  •  t*  3  * 

•  .  . 

NSMX  •  NX  -  t  •  MONO 
NSWY  .MV  -  *  *  NS  TO 
NBP  4  NX WO  *  S  ♦  | 

MBS  ■  NS I G  *8+1 
NWS  «  NSta  4  1 

•XfTE  (6*104)  LABEL*  NX. -NY*  NBX*  NBA*  NX*  «Y*-DX..OV«  (  1*  1 

t  «  3  )«  (  J*  (  6(1. J).  I  •  I*  3  It  J  t  tl  3  ) 

!E<  FLAB  ,NE*  o  )  SO  TO  8 


•ntTE  (6.104)  OUTAXE*  NtO 
TE(  fSNf  .EO*  O  )  60  TO  31 

•  *  #  0  0 

PEAO  (LTAXE.108)  (  (  B ( t « J) «  I  •  1*  NON  I*  J  •  NMS«  MBS  ) 

CK>  3  1  •  1*  NR* 

II  •  NBX  ♦  I  -  I 

so  3  j  *  i«  Nsta- 

Jt  •  MBS  4  1  -J 
B<  t  « J)  •  B(T1tJ1) 


CALL  «RfTe«<  B*  SB*  -NST6*  NS16* 


*  NXwO.  ME40( t *t )  ) 


00  4  f  •  It  NX*  a 
ft  •  !  *  4 

READ  (JTAXE.109)  (  (  N(J1«J).  J1  «  T»  1 1  ) •  J  •  I •  NY  ) 

OSAO  ( JT AXE  *  1 06  J  EL ABEL*  COOE*  TERN 
tSXI  •  3 

tE(  CODE  «EO*  AST  )  tSXt  ■  t 
IE I  CCOE  *E0*  SLASH  )  ISXt  •  1 
|E<  COOE  .EO*  OOL  3  tSXt  •  0 

TE(  TSXT  .EO*  3  )  SO  TO  SOT 

♦E(  t*T( T )  *NE*  0  ) 

J  CALL  XXfTEXt  M*  JOB*  0*  NY  -  T  t  0*  NX  -  T  *  MEAOOtX)) 

oo  6  r  •  i*  nshx 

00  6  J  •  I*  NSWY 
StJN  •  0*0 
00  9  tl  ■  1*  NBB 


6X003880  | 

6X003*30  ? 

6X003840  $ 

6XO03SBO  j 
6X003860  I 
6X003ST0  !| 

0XO038B0  I 
6X003840  I 
6X003300  1 

6X003310  I 
6X003380  I 
6X003330  | 

6X003340  I 
6XO033B0  | 

6X003360  | 

6X003370  I 
I6XO033B0  I 
6X003340  | 

60003400  I 
6X0034 TO  f 
6X003480  | 

6X003430  I 
6X003440  I 
6X003490  * 

60003460  I 
6X003470  | 

6X003460  J 
6X003440  } 

6X003900 
6X003910  * 

6X003980  j 
6X003930 
6X003940 
6X003990 
6X003960 
60003970 
6X003980 
6X003940  i 
6X003600  7 

6X0036 JO  l 
60003680  l 
6X003630  w 
6X003640  S 
60003690  I 
6X003660  f 
6X003670  J 
6X003680  1 

6X003640  I 
6X003700 
6X0037)0 
6X003780  j 
60003730  j 
6X003740 
60003790 


C-29 


oeAO  fLTAO8«t06)  8LA8€L*  COW 
t Ft  COW  ,N8.  SLASH  )  60  TO  *0* 

tSNt  •  1 

•8A©  (LTAOF.tOO)  eLAseL.  coot  - 
t8<  COW  *M€.  OOL  •  AND*  COW  (NF 
f8(  cow  ,eo,  AST  1  fSXl  •  * 

wao  ts.tost  -oex.  nlx.  ntv.  nly 
00  to  t  •  NT**  NLX 
tt  •  !  4  1  -  M8X 
00  to  J  •  MMY.-MLY 
Jt  •  J  4  1  -  NY- 
SNf  It. Jit  •  SMtt.J) 

NSHX  •  MLX  -  MMX  4  t 
NSMY  a  NLV  -  NTT  .4  1 
tori?)  •  form 

WOtTP  (8«t07)  MMX.  NLX.  NTT.  NLY 


T8(  lOJ(g)  .HS.  0  ) 

1  CALL  W»|T*St  SH«  )«•  0.  N*mY 


60003700 
60003770 
60003780 
60003700 
60003600 
60003610 
60003880 
60003830 
60003640 
60003880 
60003880 
60003870 
60003880 
60003870 
60003000 
60003010 
60003980 
60003930 
60003940 
60D03980 
60003980 
60003970 
60003980 
60003990 
60004000 
60004010 
60004880 
60004030 
68004040 
600 04 080 
60004080 
60004070 
60004060 
60004090 
60004100 
60004110 
60004 t 80 
60004130 
60004140 
60004160 
60004160 
60004170 
60004180 
60004190 
60004800 
60004810 
60004880 
60004830 
60004840 
6000*880 
60004880 
60004870 
60004880 

1.  0.  NSHX  -  1.  MFAOt 1.3)  »  60004890 


18  •  f  ♦  11  -  I 
08  3  31  •  1  *  MSS 

J*  »  J  ♦  Jt  -  1  •  •  - 

sun  •  sw  4  em.Jit  •  N<tt«J8i 
SMtt.J)  •  SUN 

tt  •  NS  Hit  ♦  1 
00  7  f  a  ft.  NX 

00  7  J  •  1 .  NSHY 
SMtt.J)  a  0.0 

•OfTE  tOUTAOfr.101  >  LAS8L.  NfO 
•Of T8  C0UTA08.il*)  OX.  07 

•OfTe  (0UTAP8U13)  ns  MX*  NSH7.  HOMO. -NS 16.  NX.  07  . 

•OlT8  f0OTAO8.11tt  t  t  Stt.J).  1  a  1 .  3  ) «  J  a  t *  3  ) 

00  71  fat.  MS  MX.  4 
tl  a  |  ♦  3 

•OfTe  (0UTA08.lt 0>  (  t  SHtJS.J).  Jt  •  t .  fl  )<  J  •  t.  NSHY  ) 

60  TO  11 


DO  9  1  a  1.  NX.  4 

1 1  •  1  ♦  3  •  •  •  *  *  •  • 

08 AO  fLTAO8.H0)  f  f  SMfft.J).  ft  •  tl  tl  It  J  a  tl  Nr  ) 

N8X  a  *  *  NX  4  1 

N6HV  a  MY  ♦  t  .... 

08A0  (LTAOe.tlO)  (  f  OMT(l.J).  1  a  1.  NTX  ).  J  •  1 .  NONY  ) 


•  AST  )  fO  TO  91 


C  -30 


ff^w.  •  •; 


Jec  •  4 
wot TP  (0*303) 

00  SO*  !  •  1*  100 

WP1TE  (OUTAPE*303) 

1P<  JEC  *tO*  4  1  WOt TP  (0*304) 
STOP 

WP1TP  (0*300)  JPC*  f 
SO  TO  lOOO 


6*004300 

tP(  ftixe  ,NS*  o  1  CACt.  TIPtMl  80004310 

6*0043*0 

call  psspio  •  spooasso 

6*004340 

tP(  tswt  .PO.  -*  i  00  TO  SO©  SPO043S0 

WOfTP  (0UTAPP*3001  0*004300 

SO  TO  tOO©  «*©043T0 

0*004300 

6*004300 

JPC  «  4  6*004400 

WPfTE  (0*303)  6*004410 

00  SO*  I  •  1*  100  6*0044*0 

•PIT*  (OUTAPP*30S1  6*004430 

1P(  JPC  *P0*  4  1  WRITE  (0*3041  6*004440 

STOP  6*004430 

6*004400 

WP1TP  (0*3001  JPC*  I  6*004470 

SO  TO  lOOO  6*004400 

6*0044*0 

VntTP  (6*301)  Pt-ASPL.*  COOP*  IP**  6*004300 

STOP  6*004310 

6*0043*0 

WP1TP  (6*30C|  ELA0EL*  CO OP*  »P*W  6*004330 

STOP  6*004340 

6*004330 

6*004360 

PO*PAT( A6*14t2*38X« 161  6*004370 

P0PPAT1 t*A6«T61  6*004360 

PO*PAT( 7*10*31  6*004340 

P0**AT(7t 1 01  6*004600 

P0**AT(*1m1 INPUT  PA»APPTP*S  — I tA0///0Xt*HlN*UT  *ATP|X4X|4*3h  X  16*004010 

*4/y0X1OMSPOOTH|N6  *ATP|X|4*3H  X  |4//tWt4wMUP6gW  OP  CA6S*X14«3m  X  |«P©O46*0 

•4//0X1 3 HO AT A  t*TE*VAl.3Xr4«t«3H  X  P4, 1 ////0X33HSPSCTPAC  SPOOTMIN6  PG»OO403O 

*AT»IX«  0(|*Jl//7X3ttO//(0Xt3«3P16*3ll  6*004040 

P0PPAT(3Pl©*31-  6*004030 

PO**AT( 1 1 A6«  A3«  A] « 18)  6POO400O 

P0*PAT(/iM03X3lKTMtS  |3  A  C0**UTaT1<JN  1NV0WVTN«/SX7HC0LUPN314*3H  T6»OO407O 

1014*1 OM  ANO  ROWS  T4*3H  TOt 4/6X83 POP  A  *WP-S*O0TPP©  *AT*tX*1  6*004060 

EGPPAT ( 40HOSENSE  SWITCH  6  UP  A*0  PWPSS  START  TO  CONTINUE/SOM  SENSE6POO404O 


W»fTP  (0*301) 
STOP 


ELAOPL*  COOP*  IP** 


WW1TP  (0*3O*|  PLA0PL*  COOP*  »P*W 

STOP 


1  SWITCH  0  DOWN  ANO  PPESS  ST APT  TO  TSPPfNATf*) 
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